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I. INTRODUCTION 


It is a well established fact that evaluations devoted to describe quantum chromodynamics (QCD) phase transitions 
need to be done in a nonperturbative fashion, such as numerically solving it on the lattice (LQCD) [1], which nowadays 
is considered the most reliable way to tackle the problem at vanishing densities. Unfortunately, LQCD is still plagued 
by the so-called sign problem [2] which prevents the method to be used to describe the transitions expected to take 
place at lower temperatures and higher densities. On the other hand, multi-loop perturbative results for many QCD 
physical quantities are within reach so that an appealing alternative would be to use them in conjunction with some 
resummation procedure in order to generate nonperturbative results. In this vein, different analytical techniques 
envisaged to combine the easiness of purely perturbative evaluations with nonperturbative optimization/resummation 
procedures have been proposed in the past decades [3]. Some of these methods are based on a reorganization of a 
given interacting Lagrangian, so that it becomes written in terms of an arbitrary mass parameter which, for massless 
theories, also works as an infrared regulator (as in hard thermal loop resummations [4, 5]). One of these approaches 
is the so-called screened perturbation theory (SPT) [6, 7], in which the variational parameter is described by a 
thermal mass. The SPT was originally proposed to describe the thermodynamics of massless scalar theories, but it 
has been later generalized so that the equation of state of thermal gauge-invariant theories [4], such as QCD, could 
also be obtained. This gauge invariant generalization known as hard thermal loop (resummed) perturbation theory 
(HTLpt) [8], has been already used to calculate QCD thermodynamic functions up to three loop order at finite values 
of the temperature and chemical potential [9, 10]. The SPT and HTLpt are actually conceptually similar to the 
so-called linear delta expansion (LDE) and optimized perturbation theory (OPT), developed earlier under various 
different names [11-13] mainly in the context of zero-temperature field theories. Within this technique, perturbative 
evaluations are performed using propagators written in terms of an arbitrary mass parameter, so that optimized 
nonperturbative results can be generated by requiring the mass parameter to satisfy a variational criterion. The two 
major problems the above mentioned methods try to solve are the poor convergence and the notoriously bad scale- 
dependence of the standard perturbative series both for the thermal mass and for the pressure at higher orders (see e.g. 
[3] for a review): not only the increasing perturbative orders show no clear sign of stability, but the scale-dependence 
worsens substantially at higher orders, at odds with what is intuitively expected for most known perturbative series 
at r = 0. Part of this bad behavior is commonly explained [3] by the unavoidable complicated interplay of soft and 
hard thermal contributions. Actually, the dynamical generation of a thermal screening mass mu ~ '/XT influences 
the relevant expansion of physical quantities, such as the pressure, which are then expressed in powers of '/X rather 
than A. In this case the predictions are, a priori, less convergent than for the T = 0 case, yet most of the interesting 
thermal physics happens at rather moderate coupling values, so that one could expect a better behavior. Despite 
the unavoidable y/X “nonperturbative” dependence, both SPT [6, 7] and OPT [14] applications to hot scalar theory 
show how these methods indeed improve the stability of the predictions when higher orders in the loop expansion are 
considered. Given the inherent technical difficulties associated with the (three loop) evaluation of the QCD pressure 
for the case of hot and dense quark matter, the recent results in [10] represent an impressive achievement. However, 
the same results exhibit a substantial increasing scale-dependence at increasing two- and three-loop orders, even for 
moderately large coupling values, which remains a surprising issue. This is more pronounced for HTLpt in QCD 
applications at the three loop level [9, 10]. While the latter are sometimes remarkably close to lattice results for 
temperatures down to T > 2Tc, for the central renormalization scale choice /i ~ 2'kT in the MS-scheme, it appears 
puzzling that a moderate scale variation of a factor 2 dramatically affects the pressure and related thermodynamical 
quantities by relative variations of order 1 or more. It is argued [10] that resumming the logarithmic dependence of 
HTLpt results may improve this problem, but as we shall explain below the missing RG invariance properties is more 
basic within the SPT/HTLpt approach. 

Recently, the standard OPT procedure at zero temperature has been modified to incorporate consistently pertur¬ 
bative renormalization group (RG) properties. It was shown to considerably improve the convergence of OPT, as 
tested for the Gross-Neveu (GN) model mass gap [15], and further used to determine with a good accuracy the basic 
QCD scale in the MS scheme, A|^, and corresponding value of the strong coupling ds[16, 17]. Very recently the 
same method has been used to estimate [18] the QCD chiral quark condensate. Here, we extend the construction to 
the case of a scalar theory with quartic interaction to show how this RG improved OPT (RGOPT) can easily cope 
with the introduction of control parameters such as the temperature, developing in detail the construction presented 
very recently in [19]. Using this textbook example we aim to illustrate how the RGOPT takes care of the scale 
dependence problems of thermal theories, those being more clearly visible within the SPT [6, 7] and HTLpt higher 
loop results [8-10], but also present in other resummation approaches. 

In a nutshell, leaving aside technical details, our basic observation is that the arbitrary variational mass introduced 
in the SPT/OPT context can (and should) be treated as any “proper” mass, from the RG viewpoint. In particular, 
this means that it involves its own standard anomalous dimension, which is unrelated to the coupling f3 function, and 
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should be incorporated consistently within RG properties. This brings crucial consequences already on the form of 
an explicitly RG-invariant pressure, prior to any further improved/resummed perturbation approaches. It also more 
clearly separates the hard and soft modes from a RG/scale dependence perspective, at least in an intermediate stage 
of the calculations. In the more standard approach both contributions are mixed up since all relevant quantities 
are expressed as a function of the coupling and temperature only. That being accepted, one realizes that most of 
the observed worsening scale-dependence at higher orders is actually due to a manifest failure of RG invariance 
for a massive theory. This is at least particularly transparent in the MS scheme largely used within SPT/HTLpt 
methods, where we show that minimally subtracting the vacuum energy divergences, without finite vacuum energy 
contributions, explicitly misses RG invariance. 

In addition, even when starting from an explicitly RG invariant perturbative expression, RG invariance is generally 
lost as a consequence of the standard (linear) modification of perturbative expansion implied in OPT/SPT, a fact 
that has been seldom appreciated in the relevant literature so far. However, RG properties can easily be restored by 
a consistent use of renormalization group properties for a massive theory, automatically incorporated in the RGOPT 
approach [17], whereby a drastically improved scale-dependence follows naturally. This opens up the possibility of 
rather simply exploiting many non-trivial SPT/HTLpt results performed up to three-loop order so far, by improving 
substantially their scale independence after applying appropriate RGOPT adaptations. 

Incidentally, other resummation approaches, like typically the so-called “two-loop-<i>-derivable” approach, (related to 
the two-particle-irreducible (2PI) method) [20, 21] are manifestly scale invariant when applied to the scalar model 
in a certain approximation. In fact one observes interesting analogies (but also important differences) between 
the first non-trival (two-loop) 2PI order and our one-loop RGOPT results, as will be elaborated in more detail 
below. However, as it happens the scale-invariance of the above-mentioned 2PI results [20, 21] is essentially due 
to a renormalization that appears peculiar to the scalar model, inspired by the properties of the 0{N) symmetric 
theory in the large-limit [22] (for which indeed the two-loop-$-derivable 2PI results become exact). Thus, to 
the best of our understanding, this scale invariance appears somewhat accidental and difficult to translate to higher 
orders and to QGD. Incidentally, the 2PI approach has been pushed even to three-loop order for the theory in 
Ref. [23], with remarkably stable results with respect to two-loop order, but the method becomes more involved, and 
a definite scale-dependence reappears, although much more moderate than in the three-loop order SPT and HTLpt 
cases. Other approaches like the nonperturbative renormalization group (NPRG) [24], should be RG invariant by 
construction. But solving the relevant NPRG equations for thermal QCD beyond approximative truncation schemes 
becomes very involved in practice. We also remark that our approach shares some qualitative features with the 
ideas invoked and the framework developed very recently in Ref. [25], in which the authors explored some form of 
massive renormalization scheme and its consequences. However, the RGOPT differs substantially in its approach as 
it incorporates by construction [17, 19] all relevant RG properties systematically order by order, relying basically on 
standard perturbation theory. Being largely based on (but not limited to) MS-scheme results, the RGOPT can easily 
be extended to any higher order calculations performed in the framework of different models, when already available, 
including thermal QGD. Therefore, it appears to us that the method presented here is conceptually simpler than 
other more sophisticated resummation approaches mentioned above. 

The paper is organized as follows. In the next section we quickly review some basic results of thermal scalar theory 
for the free energy at the relevant two-loop level. Then, in section 3 we address in some details the RG (non) invariance 
issue in the massive case, paying special attention to the MS scheme largely used in the SPT (or the similar HTLpt) 
method. In the same section we also explain how to restore the perturbative RG invariance at arbitrary orders in a 
simple fashion. Next, the resummation by optimized perturbation (OPT), with the crucial modification to maintain 
its perturbative RG invariance, RGOPT, is discussed in rather general terms in sections 4 and 5 respectively. The 
method is then illustrated in details by evaluating the free energy of a hot scalar field theory at one- and two-loop 
level in sections 6 and 7. We emphasize that all the construction developed in Secs. 3 to 7 for the (f)'^ model (for 
which we also briefly consider the large-iV case) is actually more general, most of it being straightforwardly applicable 
to thermal QCD. We occasionally mention some properties anticipated to be similar (or eventually different) in the 
QCD case. Finally, onr conclusions are presented in section 8, and one appendix deals with some technical details on 
the RG-invariant construction of counterterms. 


II. FINITE TEMPERATURE SCALAR FIELD THEORY 

We consider, as a starting point, the massive scalar field theory described by the Lagrangian density 
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where we introduce a (yet unspecified) generic mass term (m) which can be thought as a thermal mass generated 
by higher perturbative orders in an originally massless theory. One may then evaluate the free energy using known 
results from ordinary perturbation theory for the massive case. Then, up to the two-loop level the basic expression 
of the (bare) free energy is formally [3, 7, 26]: 


.Fo-^ln(p2+m2) + ^^-^^) (2.2) 

where the temperature is introduced via Matsubara’s imaginary time formalism {p^ = uj'^ + with the bosonic 
Matsubara a;„ = 2TrnT) and we have also defined 
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The divergent integrals are regulated using dimensional regularization techniques with D = 4—2e while renormalization 
is carried out in the MS-scheme. The term represents all the relevant counterterms contributions to 0(A) (see 
Ref. [7] for details). After all the mass, coupling, and vacuum energy counterterms have been consistently introduced 
to cancel the original divergences, one obtains the (MS-scheme) renormalized free energy [3, 7]: 
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where we have explicitly separated the thermal and non-thermal contributions for later convenience. Here and in all 
related renormalized expressions below, p represents the arbitrary renormalization scale introduced by dimensional 
regularization in the MS-scheme, and A = A(/i). Note carefully in Eq.(2.4) that Sq represents a possible finite vacuum 
energy term which is usually ignored, i.e. minimally set to zero in the (thermal) literature [7]. However, within our 
approach this quantity is necessarily non-zero and plays a crucial role as will become clear in the sequel. 

The standard (dimensionless) thermal integrals appearing in Eq.(2.4) are given by 


J„(x) = 4 
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where t = p/T and x = m/T. Different integrals can be easily related by employing derivatives such as 
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Also, a high-T expansion such as 
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is often useful since it represents a rather good approximation as long as a: < 1, i.e., T larger than to. 

Finally, for later comparison we recall that the thermal screening mass {mo), defined [27] by the pole of the (static) 
propagator, is obtained for weak coupling for the massless theory as a perturbative series which to lowest orders 
reads [27]: 
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III. RG INVARIANT FREE ENERGY IN MASSIVE RENORMALIZATION SCHEMES 

We now discuss the lack of RG invariance when Sq is minimally set to zero in Eq.(2.4). Remark first that to obtain 
(2.4), calculations have been performed with an arbitrary mass to in the dressed propagators (mainly to subsequently 
treat it variationally in the OPT/SPT approach), with no prejudice at this stage that it should be a thermal mass of 
order to^ ^ AT^ in the actually massless theory. Thus from the RG viewpoint, everything in (2.4) should behave as 
if it was a genuine massive theory, in particular the mass should have its standard anomalous dimension. Recall that 
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the (homogeneous) RG operator acting on a physical quantity with mass dependence, such as the free energy in the 
present case, is defined as 




d d d 
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dX 


- lm{X) 


dm 


(3.1) 


where our normalization for the fd function is 


/3(A) = 


dX 

dln/r 
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(3.2) 


while for the anomalous mass dimension it is given by 


7m (X) — 


dlnm 


dln/ii 
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(3.3) 


with [28] 


(47r)2&o = 3; (47r)^7o = (47r)% = -y , (47r)47i = -y . (3.4) 

It is easy to see that the renormalized expression, Eq. (2.4), requires a finite £o contribution to be RG-invariant, 
as one can readily see by considering the one-loop term which has an explicit In/i dependence. Thus, acting with 
the RG operator, Eq. (3.1), on the RHS of Eq. (2.4) gives a non-zero contribution of order 0(1): —(l/2)m"^, which 
is not compensated by terms in Eq. (3.1) coming from the lowest orders in /3(X)d\ or m7m(X)dm oc Xm'^, those 
being at least of next order 0(A). This is a manifestation of the fact that (perturbative) RG invariance generally 
occurs from cancellations between terms coming from RG coefficients at order and the explicit fx dependence at 
the next order This can also be understood alternatively by considering solely the original bare contribution 

to the free energy: although the latter only depends on the RG-invariant bare mass and coupling too,Ao (and on 
2e = 4 — H in dimensional regularization), its finite part is not a priori separately RG invariant. In other words 
for a massive theory the T-independent vacuum energy divergences cannot be absorbed by an arbitrary redefinition 
of the vacuum energy without spoiling RG-invariance. Now, as we recall below the subtraction needed to recover 
RG invariance is perturbatively well-defined and easy to construct order by order. The vacuum energy gets its own 
anomalous dimension which, within dimensional regularization, is essentially determined by the coefficients of the 
poles in 2e = 4 — H, stemming from the remaining divergences once the mass and the coupling have been properly 
renormalized. This procedure had been exploited in an earlier application of the OPT to evaluate the vacuum energy 
of the Gross-Neveu (massive) model [29] and then extended to the QCD case [30, 31]. Similarly, a well-known related 
result is that the Coleman-Weinberg effective potential for a general massive theory is not RG invariant without finite 
“vacuum energy” terms independent of the fields, as was originally carried out in Ref. [32] and in the MS-scheme in 
the context of RG-improvements of the effective potential [33]. Indeed, for the 0{N) model, the vacuum energy 
anomalous dimension has even been computed up to four- and five-loop order in Ref. [34]. 

However, to the best of our knowledge, this point appears to have been overlooked in the context of thermal theories. 
In applications of improved/resummed massive perturbation schemes such as SPT [6], HTLpt [9], and the standard 
OPT [14] the calculations are mostly performed within the MS-scheme and the T = 0 vacuum energy divergence 
is minimally cancelled out by appropriate (zero point) counterterms but missing out those extra finite subtractions 
required by RG properties. In fact, as far as the purely perturbative massless theory is concerned, the only mass 
is actually a thermal mass: ~ AT^, so that the lack of RG invariance pointed out above is rather postponed 

to higher (three-loop) perturbative order A^, where it plainly resurfaces. Within the SPT, or the similar HTLpt, 
approaches the variational mass parameter is similarly perturbatively of order wf ~ XT^. It is thus not surprising 
that the scale dependence observed within SPT/HTLpt results appears to worsen at higher orders [7, 9]. But more 
generally one wishes to use the nonperturbative mass gap resulting from such variational approaches possibly beyond 
standard perturbation for moderately large coupling values, as can be relevant near a critical temperature. Thus, the 
lack of RG invariance appears more serious since as we recall in next section (see Eq. (4.1), the variational procedure 
formally treats the mass to be of the same perturbative order as the lowest order considered contributions, like e.g. 
the “hard” thermal one-loop contributions of order ~ A° T'^. Moreover, in the standard procedure one makes the 
arbitrary renormalization scale pL effectively temperature dependent by choosing p, ~ 2ttT, such as to avoid large 
ln/x/(27rT) contributions coming from the remnant scale-dependence. In this way the pressure can be studied as a 
function of T/Tc in QCD applications, where Tc is related to the basic QCD Aqcd, e.g. in the MS-scheme. But if 
the scale dependence appears not much reliable at higher loop orders, one may question as well the reliability of the 
corresponding T/Tc dependence of the pressure, even for the well-motivated central p = 2'kT prescription. 



6 


While those issues in MS or related schemes may perhaps not explain at once all the problems that thermal theories 
face with perturbative expansions at increasing coupling, a part of those problems are likely to be reduced if one adopts 
from the beginning a prescription fully consistent with RG properties. This problem appears partly circumvented 
(but are actually rather delayed to higher orders) in thermal perturbative calculations performed in some other renor¬ 
malization schemes, where the zero-point energy, Tq{T = 0), is subtracted for convenience prior to any subsequent 
calculations. Indeed, subtracting the T = 0 contribution from Eq. (2.4) washes out all the first /i-dependent terms, 
making scale-independence (trivially) satisfied at one-loop order, and (less trivially) at the two loop level as well. But 
then the one-loop result becomes also trivial, with the only left contribution being the pure thermal, third term in 
the RHS of Eq.(2.4). So, there are no possible optimized solutions of the OPT/SPT/HTLpt form, which can only be 
obtained at the two-loop level. Moreover, applying subsequently the standard SPT/OPT procedure anyway spoils 
RG invariance. In any case it is most convenient to have a prescription generically valid both for zero and finite 
temperatures, so that subtracting T = 0 contributions is not satisfactory for a more general framework. Accordingly, 
the subtraction procedure we will consider next only depends on T = 0 contributions but is generically valid also 
for T ^ 0. Moreover, a remarkable consequence is that the subsequent mass optimization, as implied by RGOPT, 
will give a nontrivial solution already at the one-loop order, and very similar to what is normally obtained at two- 
loop order with the other mentioned resummation schemes (SPT/OPT, HTLPT, 2PI,...), as we will examine in detail. 


Following Refs. [16, 17, 29-31] the easiest way to construct an RG-invariant finite vacuum energy is to determine 
So order by order as a perturbative series from the reminder of acting with Eq. (3.1) on the non RG-invariant finite 
part of Eq. (2.4): 

= -Remnant(A,m) = - fj,-^[J^o{So = 0)|fl„ite] (3.5) 

d^j, dfi 

where the RHS of (3.5) thus defines the anomalous dimension of the vacuum energy. As above mentioned, it is easy 
to see that, as a perturbative series, Sq has the convenient form in MS or similar schemes 

^o(A,to) =VsfeA'" , (3.6) 

fc >0 

where the constant coefficients Sk are perturbatively determined order by order, being essentially determined by 
the coefficients of the (single) powers of In /i term at order fc -I- 1 (or equivalently by the single poles in 1/e of the 
unrenormalized expression) [17, 30]. This procedure leaves non RG-invariant remnant terms of perturbative higher 
orders to be cared for similarly once higher order terms are considered. The apparently odd divergent behavior for 
A —>■ 0 of this first order term is actually not a problem since, as we will see explicitly, it completely disappears from 
the final results. 

We stress that all the previous considerations, being only dependent on the renormalization procedure, do not de¬ 
pend on the thermal contribution so that, at arbitrary perturbative orders, the subtraction function represented by 
Eq.(3.6) can be determined simply from the T = 0 contributions only Note that Eq. (3.6) is not the only possible 
subtraction form in general, but a very convenient one in MS or related schemes to proceed systematically at higher 
orders. Moreover, it is particularly convenient once introducing below the RGOPT modification of perturbation, 
since the 1/A term will be responsible for a non trivial RGOPT solution already at one-loop order. 


Explicitly, one finds the RG-invariant form of Eq. (2.4) up to two-loop order with a non-trivial Sq given by 

So = -m* + Si + S 2 A + 0(A") 

After some algebra one obtains 

*0 = 2(60-470) ^ 

9677^ (6o-1287r=^((l-14si)7i-so(62-472))-41 23-136C[3] ^ , n onn 

*2— 12288 7i-"‘(6o-I-47o) “ 48n-7r2 — O.UiOUU , 


(3.7) 


_ (61-471) _ , 

^ 870(^0—470) ’ 


(3.8) 


where the explicit RG dependence in the intermediate terms emphasizes the more general form of these results, while 
the last terms are specific to the N = 1 (j)'^ theory. To derive S 2 according to the previous discussion we had to use 


^ However, when the T = 0 and the T 7 ^ 0 contributions are not explicitly separated, like in the case of two- and three-loop HTLpt [9], 
due to the systematic m/T expansion making such involved calculations tractable, caution will be needed to expand at a sufficient order 
in m/T so as to get all the relevant terms of the same perturbative order to construct the corresponding subtractions in Eq. (3.6). 
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the (T = 0) In^ coefficient at three-loop order given e.g. in Ref. [7]. 

One may equivalently derive the finite subtraction in Eq. (3.7) in an alternative manner by RG invariance consider¬ 
ations solely on the bare expression of the free energy. For completeness, this is presented in the appendix. Instead 
of minimally subtracting the bare vacuum energy divergence, an RG-invariant counterterm can be added to cancel 
the remnant divergences, and necessarily incorporates also the same finite subtraction terms in Eq. (3.8). As a 
non-trivial crosscheck of our calculation, let us note that now acting with the RG operator, Eq. (3.5), on the results 
given by Eqs. (3.7) and (3.8) one recovers, for = 1, the results up to of the anomalous dimension /3«(A) (with 
fidSo/d^j, = 2m^/3„(A)) which has been calculated up to four and five loops for arbitrary N in Ref. [34]. (Actually 
we could have used directly the results in [34] in the present scalar model case to derive the Sk in Eq. (3.8), but the 
above derivation using the basic available perturbative expressions shows precisely how to proceed for an arbitrary 
theory, where the vacuum energy anomalous dimension may not always be explicitly available.) 

One sees that the subtraction with Sk, explicitly depending on RG coefficients, incorporates a non-trivial RG- 
dependence already at first (one-loop) order, only depending on already known one-loop standard RG coefficients. 
This result has important consequences for the subsequent OPT application. Now, there is a subtlety at this stage: 
while the Sk subtraction terms are strictly necessary to recover RG invariance at order A^, i.e. up to neglected 
terms, they enter the free energy expression at order as Eq. (3.6) indicates. For instance only sq is needed to 
recover RG-invariance at one-loop 0(1), but the next term si is of order 0(1), so strictly speaking si should also be 
included in the full “one-loop” free energy results. This appears as a complication a priori, meaning that at order k 
one needs in principle the more demanding information from (the In/i coefficient of) perturbative order fc-l- 1. On the 
other hand since RG invariance is constructed perturbatively, one may expect that the simplest minimal prescription 
of keeping only the Sk terms at order A^ should already be a good enough approximation. Accordingly, we mainly 
consider below the simplest prescription but also examine both prescriptions, indicating the differences whenever 
relevant. We will see that, after the modification of perturbative expansion implied by RGOPT, incorporating the 
higher order Sfc+i at order A^ makes no crucial differences, even at one-loop order, the resummation results being 
not very sensitive to such purely perturbative variations. The same stability with respect to such variations was 
also observed at vanishing temperature in Ref. [17] (where those different prescriptions were incorporated within the 
intrinsical theoretical errors of the method). 

It should be clear from the previous derivation that by construction the subtraction terms make the free energy 
perturbatively RG-invariant. But just to crosscheck it in a more pedestrian way, let us now reexamine the result at 
one-loop order, with the —SQrn^/X subtraction included in Eq. (2.4), using the standard one-loop RG running coupling 
and mass. These are given from integrating respectively Eq. (3.2), (3.3) which yields the standard textbook result: 

A(^) = A(^o) - feoA(Mo)ln^^ ~ A(^o) + Y^^(mo) In ^-f 0(A^)^ , (3.9) 

to(^) = m(/ro) (^1 - 6 oA(/ro) In~ m(/ro) (^1-f-^A(/ro) In-H e>(A^)^ , (3.10) 

V doj \ 327r^ ^0 / 

where A(/io) and m(/xo) are the coupling and mass at some arbitrary reference scale fio. Expanding Eq.(2.4) to first 
order in A (i.e. to order A°) one obtains: 

- 3 ^ - i(41n/i-t3H- ) + - j^\n^i) + 0{X) -t thermal part 

^niHdo)[-^-l+0{X)] , (3.11) 

(where • • • stands for /x-independent terms). This relation explicitly displays the cancellation, up to terms of higher 
order 0{X), of all In/x contributions coming respectively from the original one in Eq. (2.4), and from the running 
coupling and mass (third and fourth terms respectively in Eq.(3.11)). We have neglected above the thermal con¬ 
tributions, but including these does not alter the results since the exact T-dependent contribution does not depend 
explicitly on /x. Alternatively, when the high-T expansion is considered, the explicit In/x/m in Eq. (2.4) is replaced 
by a ln/x/( 47 rT) with the same coefficient consistently. It is instructive to push this exercise a step further now 
incorporating in the free-energy (2.4), restricted at one-loop order, all thermal contributions in the high-T limit, from 
Eq.(2.7), and plugging into the resulting expression m = mo, the standard perturbative thermal mass, Eq.(2.8) also 
restricted at first order. Then, the expression only depends on the coupling, and using Eq. (3.9) gives the result: 

+ 0 (AMn^) 


(4^)2j-o ~ T^ 


(3.12) 



where all fi dependence has been cancelled out up to order A^, since ^ A^. This was expected since the subtraction 
takes care of the T = 0 lowest orders /x-dependence as shown previously, while the thermal contribution oc jQ{m/T) in 
Eq. (2.4) does not depend on the scale. More interestingly, one recognizes from Eq. (3.12) the standard perturbative 
expansion for the pressure which, upon using the common normalization with the free gas pressure Pq = tt^T'^/QO 
and Ao = 4!g^, can be expressed as: 


_P ^ 15^ 

Pq 8 TT^ 2 TT^ 


(3.13) 


In particular, the term —Sq/A is incorporated in the final result Eq. (3.12) once using Eq.(2.8) since —Somjj/X ~ 0(A). 
The latter gives a contribution (15/8) to P/Pq, making this complete one-loop expression consistent with the stan¬ 
dard perturbative expression of the pressure up to order A^^^, while the original (unsubtracted) perturbative one-loop 
expression is not, giving an expansion similar to (3.13) but with a twice too large term —(15/4) g^. However, this 
agreement is merely an accident of one-loop order: at two-loop order, including the next order subtraction term from 
(3.7) does not give the correct massless perturbative pressure when replacing m by itid- (In particular the subtraction 
—som'^/X happens to be exactly cancelled by the hard contribution of order A, oc XT^). But this is not surprising, 
since more generally the perturbative massless pressure [27] cannot be obtained consistently from simply replacing 
the perturbative thermal mass (2.8) in the expression of the massive pressure. However, when the mass is traded 
for a variational parameter, like in the OPT/SPT resummation approaches to be recalled next, one may recover 
the massless pressure results under specific conditions. This is made possible because the OPT construction and 
mass optimization drastically modifies the massive contributions as compared with the original perturbative expansion. 


To summarize this section, starting with Eq. (2.4) including £q from Eq. (3.7) to obtain a perturbatively RG- 
invariant free energy provides a sound basis for more elaborate resummation procedures like the OPT method to be 
addressed next. As above mentioned, an extra advantage of the subtraction terms (3.7) starting with 1/A, is that the 
optimization procedure will provide a non-trivial mass m(A) already at the lowest one-loop order, in contrast with the 
standard SPT and HTLpt approaches (where at one-loop order the mass optimization gives a trivial A-independent 
solution [6-8]). This is a welcome feature specially for more involved theories like thermal QCD where higher order 
contributions are challenging to evaluate, and a comparison between successive orders is crucial to establish the 
stability of the resummation results. 


IV. OPTIMIZED PERTURBATION THEORY (OPT) 

The basic feature of the optimized perturbation theory (OPT) (appearing also under different names and varia¬ 
tions [11-13]), is to introduce an extra parameter 0 < 5 < 1, which interpolates between Cfree and Cmt in Eq.(2.1), so 
that the mass m is traded for an arbitrary trial parameter. This is perturbatively equivalent to taking any standard 
perturbative expansions in A(/x), after renormalization, reexpanded in powers of S after substituting: 

m —>■ m (1 — 5)“, A —>■ (5 A . (4.1) 

This procedure is consistent with renormalizability [29, 35, 36] and gauge invariance [31], whenever the latter is 
relevant, provided of course that the above redefinition of the coupling is performed consistently for all interaction 
terms and counterterms appropriate for renormalizability and gauge invariance in a given theory^. Note that in 
Eq. (4.1) we have introduced an extra parameter, a, to reflect a priori a certain freedom in the interpolation form. 
As will be demonstrated below this parameter plays an essential role within our method for allowing compelling 
constraints to be imposed. Applying Eq. (4.1) to some given renormalized perturbative expansion for a physical 
quantity, P(m, A), reexpanding in d to order k, and taking afterwards the <5 —>■ 1 limit (to recover the original 
massless theory) leaves a remnant m-dependence at any finite J^-order. The arbitrary mass parameter m is then 
most conveniently fixed by a variational optimization prescription known as the principle of minimal sensitivity [13] 

^pW(m,A,5=l)l™^A = 0, (4.2) 

am 


^ Contrary to what is sometimes claimed (and worked out) in the literature, the OPT/SPT/HTLpt does not need any extra counterterms 
besides the standard ones of the corresponding massive theory: in particular all seemingly new divergences generated at arbitrary orders 
from the first replacement in Eq. (4.1) are evidently related to the single standard mass counterterm. Moreover at arbitrary orders, 
temperature-dependent divergences and associated counterterms should not appear, as expected from general principles, provided that 
one keeps the mass as an arbitrary parameter and carefully subtract all nested subdivergences until renormalization has been completed, 
before using a gap-equation giving perturbatively m oc T \/A- 
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thus determining a nontrivial optimized mass rh{X), with nonperturbative A-dependence, realizing dimensional trans¬ 
mutation (more precisely, e.g. for asymptotically free theories at vanishing temperatures, the optimized mass is 
automatically of the order of the basic scale A ~ in contrast with the original vanishing mass). 

In simpler {D = 1) models, at vanishing temperatures, the procedure may be seen as a particular case of “order- 
dependent mapping” [12], which has been proven [37] to converge exponentially fast for the D = 1 (f>^ oscillator 
energy levels. For higher dimensional D > 1 renormalizable models, no rigorous convergence proof exists, although 
the OPT was shown to partially damp the factorially divergent (infrared renormalons) perturbative behavior at large 
orders [38]. Nevertheless, this technique can give rather successful approximations to some nonperturbative quan¬ 
tities beyond the large-(or mean field) approximations in a large variety of physical situations which include the 
study of phase transitions within condensed matter related renormalizable models [39-42] as well as within QCD non 
renormalizable effective models [43]. 

We emphasize that at finite temperatures the very same basic idea has been exploited by the SPT [5, 6]/HTLpt [8] 
method, where in this thermal context the screening thermal mass is treated as an arbitrary variational parameter, 
and in Eq. (4.2) P also depends on T like e.g., Eq. (2.4). 


V. RENORMALIZATION GROUP COMPATIBILITY OF OPT 

In most previous standard OPT (or similarly SPT and HTLpt) applications, the so-called linear ^-expansion is used, 
assuming a = 1/2 (i.e. rn? —>■ m^(l —5) for a scalar mass) in Eq.(4.1) mainly for simplicity and economy of parameters 
while the more recent approach, developed in Refs. [15-17], differs in two crucial aspects which turn out to drastically 
improve the convergence. First, it introduces a straightforward marriage between OPT and renormalization group 
(RG) properties, by requiring the (d-modified) expansion to satisfy, in addition to the OPT Eq.(4.2), a standard RG 
equation: 


= = (5.1) 

a fx 

where the RG operator was defined in Eq. (5.1). Moreover, once combined with Eq.(4.2), the RG equation takes the 
reduced massless form: 


P^'^\m,X,6 = 1) = 0 . (5.2) 

Therefore, Eqs. (5.2) and (4.2) if used together, completely determine optimized m = rh and g = g “variational” fixed 
point values. 

Since interaction and free terms from the original perturbative series are rather drastically reshuffled by the modifi¬ 
cation implied by Eq. (4.1), the RG invariance is in general no longer perturbatively satisfied, even when the original 
perturbative series is RG-invariant prior to performing (4.1). This spoiled RG invariance has to be restored in some 
manner, and thus Eq. (5.1) gives a nontrivial additional constraint. This feature has been seldom appreciated and 
considered in many former applications of the (5-expansion/OPT method to renormalizable theories (perhaps in part 
because in many analyses with more elaborated theories the OPT is restricted to first order, where RG improve¬ 
ments are supposed to play a minor role). This important role of RG properties was recognized much earlier in Refs. 
[29, 30] where to recover the RG consistency the standard linear (5-expansion was resummed to all orders. Indeed, 
this resummation can be done, at least for the pure RG dependence up to two-loop, but the result comes as a rather 
involved integral representation, not practically intuitive and making difficult to perform the mass optimization or to 
generalize to other physical quantities and other models. In contrast, the purely perturbative procedure together with 
Eq.(5.1) appears a considerable shortcut, straightforward to apply to any model. Intuitively, just as the stationary 
point solutions from Eq. (4.2) are expected to give sensible approximations, at successive orders, to the actually 
massless theory, one similarly expects that combining the latter with the RG solutions should further give a sensible 
sequence of best approximations to the exactly scale invariant all order results. 

Still, a well-known drawback of the standard OPT approach is that, beyond lowest order, solving Eq. (4.2) generally 
gives more and more solutions at increasing orders, some of which are likely to be complex-valued. More generally, 
without some insight on the nonperturbative behavior of the solutions, it may be difficult to select the right one, 
and the unphysical complex-valued optimized solutions at higher orders are embarrassing. This is incidentally a 
problem encountered first at three-loop order in SPT [7] and HTLpt applications to QCD [9]. The mass optimization 
is then replaced by alternative prescriptions, most often using simply the purely perturbative screening mass, but 
accordingly loosing a more nonperturbative ingredient from the optimized mass. But RG considerations also provide 
a possible way out, which is the second main difference and new feature of the present RGOPT version. For QCD a 


d d 
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compelling selection criterion was proposed, in which only the solution(s) continuously matching the standard pertur¬ 
bative (asymptotically free) RG behavior for vanishing coupling are retained [16, 17]. This prescription can easily be 
generalized to any model, like the non asymptotically free (AF) (j)'^ theory, by similarly requiring to asymptotically 
match the solutions to the standard perturbative behavior for small coupling, namely for fixed m and arbitrary scale 
M: 


A(Ai < m) ~ In + O ^ . (5.3) 

At zero temperature this turns out to give a unique solution for both the RG and OPT equations, up to rather high 
orders. An additional welcome feature is that by requiring at least one RG solution to fulfill Eq. (5.3) leads to a 
strong necessary condition on the basic interpolation, Eq. (4.1), uniquely determining a from the universal (scheme- 
independent) first order RG coefficients: a = ■jo/bo, as we derive in more detail below. A connection of the OPT 
exponent a with RG anomalous dimensions/critical exponents had also been established in a very different context, 
in the D = Z (fA model for the Bose-Einstein condensate (BEG) critical temperature shift by two independent OPT 
approaches [41, 42], where it also led to real OPT solutions [42]. However, AF-compatibility and reality of solutions 
can appear to be mutually exclusive beyond lowest order, depending on the particular model. A simple way out is 
to further exploit the RG freedom, considering a perturbative renormalization scheme change to attempt to recover 
RGOPT solutions both AF-compatible and real [17]. We will see that this extra complication is not even necessary in 
the present case, where (at least up to the two-loop order) the RG-compatible solutions remain real for a large range 
of relevant values of the coupling and temperature. All these features are easy to generalize at finite temperatures 
due to the fact that RG properties are essentially determined by the divergence structure of the T = 0 part. So, the 
only complication is technical since at finite temperature the previous Eqs. (4.2), (5.1), and (5.2) come with an extra 
T dependence. Let us now illustrate explicitly all those features by evaluating the RGOPT modification of the free 
energy of a thermal scalar field. 


VI. 1-LOOP, 0{5°) 
A. T = 0 


Let us first truncate Eq. (2.4) at strict one-loop order, and first restricting to T = 0 which is simpler and sufficient 
to determine the RG-exponent a in Eq.(4.1). We have 


( 47 r )2 = 0) = ( 3 -f 2 In A_ 


A 




( 6 . 1 ) 


where the superscript RGI emphasizes the (perturbative) RG invariance of this quantity. At this order the calculation 
is elementary so it can best illustrate the main steps. Now, applying Eq. (4.1), performing the i5-expansion to order 
(5° consistently, and taking afterwards (5 —>■ 1, gives: 


(47r)2j*Gi(,5°,(5 = 1) = 771^ 


-y(l -4a) - i (3-f 21n-^ 
A 8 V 


( 6 . 2 ) 


Note that the 0{1) term remained unmodified (this is a general property of OPT: expanding at order 5^ and taking 
(5 —>■ 1 leaves the order A* term unaffected due to the screening from A —>■ 6X). Then, requiring Eq. (6.2) to be 
perturbatively RG-invariant after this modification of the perturbative series, i.e. applying the RG Eq. (5.2), on gets 


1- — a +0(A) 
70 


= 0 , 


which uniquely fixes 


(6.3) 


_ To _ 1 
“ 5o “ 6 ’ 


(6.4) 


where, in Eqs. (6.3) and (6.4), we made the RG coefficient dependence explicit to emphasize the generality of these 
results. At this point several remarks are in order: 
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• 1) the very same result, Eq. (6.4), was obtained [16, 17] (up to a trivial bo difference of normalization by a factor 
2), while considering the RGOPT for QCD (with appropriate QCD values for those RG coefficients). The expo¬ 
nent a is universal for a given model, in the sense that it only depends on the first-order RG coefficients, which 
are renormalization scheme independent. Furthermore, at vanishing temperature, Eq. (6.4) greatly improves 
the convergence of the procedure at higher orders: considering only the Hrst RG coefficients 6o:7o dependence 
(i.e., neglecting higher RG orders and non-RG terms), it gives the exact nonperturbatively resummed result at 
the first 6 order and any successive order [17]. This is not the case for a = 1/2 (for a scalar mass), where the 
convergence appears very slow, if any. 

• 2) The standard linear J-expansion interpolation, widely used for zero temperature models and for SPT/HTLpt, 
takes a = 1/2 [7, 11, 14, 36]: 

—>■ m^(l — 5) , (6.5) 

thus our RG-compatible exponent (6.4) is three times smaller Indeed the standard OPT/SPT interpolation 
(6.5) leads to an unmatched RG equation, while the OPT equation (4.2) solved for A(m) (or equivalently for 
m(A) gives, in the perturbative regime /i <C m: 

^ , , 167r^ 

A(/r<m)-’ (6-6) 

in clear contradiction with the true running by a wrong overall sign plus a factor three too large. 

• 3) At the very first non-trivial order, once having fixed a = jo/bo the RG equation is satisfied and thus does 
not give further constraint. We will see that at the next and higher orders in 6, Eq. (6.4) is always required for 
the RG equation to have at least one solution matching Eq. (5.3). In addition, it also fixes A in terms of the 
other parameters (m, and the only remaining parameter /i/T when considering the thermal part). 

Let us next consider the other OPT constraint given by Eq. (4.2). Still neglecting the thermal part, and pulling out 
an overall factor yields 


boX 




= 0 . 


(6.7) 


One readily remarks the explicit exact scale-invariance of Eq.(6.7), thus of its solution, provided that one uses for 
A = A(/x) the exact (one-loop resummed) running in Eq. (3.9), since the expression l/A(/i) -I- 6 oln/x is explicitly 
/i-independent. Letting apart the trivial m = 0 solution, Eq.(6.7) has the more interesting solution 

, ( 6 . 8 ) 


which is seen to be compatible, for A —>■ 0, with the perturbative first order standard RG solution obtained from 
solving Eq.(3.2) for X{ii) at first order. Namely, for fixed m and arbitrarily small /i, it exhibits infrared freedom: 
A(/x <C to) ~ {bo \n{m/fij)~^. 

Moreover, plugging Eq. (6.8) within the modified vacuum energy expression (6.2) for the case of vanishing tempera¬ 
tures, gives a remarkably simple result: 


Jo (to. A) 



(6.9) 


Despite its somewhat trivial look, Eq. (6.9) represents a nonperturbative result, in the sense that it only involves 
the expression of to^ from Eq. (6.8). Apart from the one in Eq. (6.8) all A-dependence disappeared (in particular the 
—So/A term has consistently disappeared upon using the OPT gap equation (6.7)). Therefore, Eq. (6.9) gives a non¬ 
trivial T = 0 (negative) vacuum energy contribution and resembles much the large-iV nonperturbative result, up to 
appropriate identification of boX, but here obtained from the OPT. When higher RG orders and non-RG contributions 
are included, they spoil this simple form result [17], as does also the thermal part, that we will consider next. In 
these cases there are remnant coupling dependences, once having used the solution of Eq. (4.2) within the physical 
expression of the vacuum energy. 


® In QCD e.g for 3 light quark flavors a = 7o/(2f)o) = 4/9 is also substantially smaller than the linear case value a = 1 for fermion masses. 
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B. T / 0 


Let us now consider the thermal contributions in Eq. (2.4) still at one-loop order. After performing the (5-expansion 
to the lowest corresponding order we obtain the T ^ 0 free energy, similarly as in (6.2): 


^ 0 , 5 °, b=\) = m^ 


1 

2bo A 


1 

8 


3 -h 2 In ■ 



( 6 . 10 ) 


where the first term —som^(l —4a)/A = —m‘^/{2bo A) is the only one affected by (4.1) at this lowest order. Calculations 
are slightly more involved than for T = 0 but note that the (5-expansion and subsequent OPT minimization equation 
involve successive derivatives of the thermal function J„(m/T). For more generality and to make contact with various 
other resummation methods, it turns out to be particularly convenient to express all our RGOPT results in terms of 
the one-loop renormalized self-energy, including all thermal dependence: 



1 

p2 _|_ ^2 




7oA 





( 6 . 11 ) 


(where for completeness the mass counterterm reads reads /{2e)). This simple factorization is possible for 

the scalar (j)^ model up to the two-loop level, because the two-loop contribution (the last order A term in Eq. (2.4)) 
factorizes as the square of one-loop expressions (i.e. graphs with a different “nested” topology only appear at the 
three-loop level for the \(j>^ interactions). 

Then, noting that f^ln(p^ + iv?) = 2E/j/A, the exact solution of the OPT Eq. (4.2) can easily be written in 
the form of a self-consistent “gap” equation for fh: 


(47r)^ boT,R = b 


°2 



+ T‘^Ji 



( 6 . 12 ) 


which like the T = 0 previous case, is exactly scale-invariant by construction, as we will illustrate more explicitly 
below. 


1. Digression: connection with large-N and 2PI results 

As an important digression, we point out that Eq. (6.12) is recognized as the very same solution obtained for the 
large-A^ 0{N) (f)^ model in Ref. [22], upon appropriate 5o definition for the large-A^ case. Indeed, in the leading 1/A^ 
approximation, the only contributing graphs have the one-loop structure, and the mass-gap equation can be solved 
exactly. More precisely it is easily found from the arbitrary N RG coefficients given e.g. in Ref. [28] that in our 
normalization. 


(47r)2 bo{N) = ^ , (47r)^ 7o(Af) = 

so that in particular we have for the crucial exponent in (4.1): 



1 /Af + 2A 

2 \N + 8) 


—>• - for A^ —>• oo 
2 


(6.13) 


(6.14) 


for which values all our previous construction, and the corresponding mass gap equation in (6.12), reproduce exactly 
the large N results in [22]. Note in particular that a = 1/6 for A^ = 1 while accidentally the large N value of a = 1/2 
is the standard linear one, but here being fully consistent with RG properties. The fact that the one-loop RGOPT 
reproduces exactly the large N result can be seen as the finite temperature analog of similar RGOPT properties [15] 
obtained for the large N limit of the GN model. 

Similarly, Eq. (6.12) is also recognized as the very same form of mass gap solution obtained in the 2PI formalism but 
at two-loop order [21] (or also in the tadpole approximation for the self-energy [20]), except that in [21] the correct 
bo = 3/(167r^) is effectively replaced by bo/3 because, as explained by the authors, only one channel out of three 
is taken into account at this level of the 2PI approximation, similarly to the leading 1/N approximation. We will 
come back below on this apparent bo value mismatch when discussing the perturbative reexpansion of the pressure to 
make contact with standard perturbation results. In fact, the analogy with Ref. [21] goes further, in particular their 
expression of the mass gap solution is also exactly scale invariant at two-loop order, and the free energy involves, after 
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renormalization, a term —to'^/( 2A), once again identical to our subtraction terms —m'^/{32Tr^boX) in Eq. (6.2), when 
taking &o &o/3. But the origin of this 1/A term in ref. [21] is very different, emerging upon using the gap equation. 
Incidentally, the scale-invariance of the 2PI results is essentially due to a renormalization procedure that is peculiar 
to the scalar model, inspired by the large N limit [22], which thus appears to us as being accidental, and limited to 
the first non-trivial two-loop level. In contrast the mass gap (6.12) is obtained already at the one-loop level, and it 
should be clear from the previous construction that the RGOPT systematic procedure is applicable in any model at 
arbitrary orders. 


2. One-loop RGOPT mass gap and pressure solutions 


Although Eq. (6.12) may easily be solved numerically, it is instructive to consider next the approximation given by 
the high-T expansion, which is very precise as long as T > m: in fact this condition can be easily checked a posteriori 
considering the optimized solution m. At one-loop level, it turns out that the optimized mass always satisfies this 
criterion for all the relevant range of coupling values. Therefore, we will not need to solve the exact Eq. (6.12) for all 
practical purposes. In the high-T approximation (2.7), the one-loop order OPT Eq. (4.2) produces (discarding the 
trivial solution m = 0) a simple quadratic equation for m: 


bnX 


-f Lt 


■2tt X — 2— = 0 , 
o 


(6.15) 


where m = xT and we defined for shorthand notations the /i/T dependent part Lp = ln[/i/(47rT)e”’'®]. As already 
explained above the RG Eq. (5.2) reduces to Eq. (6.3) which is already satisfied for Eq. (6.4), thus it gives no additional 
constraint. 

Solving Eq. (6.15) gives two real solutions, but one is clearly unphysical, giving m < 0 for any A. The other unique 
physical solution is thus 





(6.16) 


where for more generality we kept the dependence on 6o explicit. Despite the apparently minor modification of the 
series represented by Eq. (6.2) at this first order, the solution given by Eq. (6.16) has clearly a nonperturbative 
dependence on A. We stress that the variational mass (6.16) is unrelated to the physical screening mass [27] in 
Eq. (2.8), and thus has no reason to reproduce the latter. Moreover, as anticipated, Eq. (6.16) is strictly exactly 
scale-invariant, provided that one uses for A = A(/i) the exact (one-loop resummed) running in Eq. (3.9) (now with 
/i being temperature-dependent as usual), since the expression l/A(/r) -|- boLp = l/A(/i) -I- 6o In/r -I- • • • is explicitly 
/i-independent. In other words the mass gap in (6.16) actually only depends on the single parameter boX{^o), where 
/To is some reference scale, typically /io = 2TrT. 


Before we proceed, it is worth to comment a little more on this result: recall that prior to the d-expansion, the 
basic one-loop expression (first line in Eq.(2.4)) with the first term sq in the subtraction £q, is by construction RG- 
invariant to one-loop order, i.e. up to neglected higher order terms 0{X). But what is more remarkable is that the 
optimal mass resulting from solving (4.2) is exactly scale-invariant to all orders (of course “all orders” but neglecting 
genuine higher orders in the running coupling, i.e. keeping only the 6o dependence to all orders). This is a direct 
consequence of the value a = 'yo/bo in the interpolating relation, Eq.(4.1). This result is the finite temperature analog 
of what was similarly obtained generically for zero temperature QCD in Ref. [17]: namely, “all-order” (one-loop RG) 
resummed results are correctly obtained by the very first RGOPT 6 order. Indeed, since this is a generic result, 
we anticipate that applying the same procedure to thermal QGD will give similar one-loop results, with an OPT 
equation and solution very similar to Eqs. (6.15) and (6.16) up to obvious changes in some factors, but also exhibiting 
exact scale-invariance). However, this exact scale-invariance is due to the form of the exact one-loop running of the 
coupling, perfectly matching Eqs. (6.15), (6.16), which does not generalize once including higher RG orders in the (3 
function and non-RG dependence at higher orders. As we examine in next section, at the two-loop 0{X) order, the 
scale invariance resulting from RGOPT extends beyond the two-loop perturbative order at which it is imposed by 
construction, but a (moderate) scale dependence reappears unavoidably at a finite higher perturbative order, precisely 
at order A^, thus one order higher than naively expected. 

To proceed one may expand Eq. (6.16) perturbatively, which is easily seen to be an expansion in '/X, as expected. 
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One then finds: 


~Y' 


TT 



^(3 - 2Lt)(&oA)3/" + Lt (boXf + ■■■ 


(6.17) 


where we kept the bo dependence explicit on purpose. 

As an important side remark, np to now we have considered the simplest minimal prescription of incorporating only 
the —(so/A) subtraction at one-loop order, strictly necessary for recovering perturbative RG invariance. It is thus 
opportune to mention what is changing if incorporating the next order subtraction si 0 term from (3.7), being 
formally also of one-loop order. In fact this simply amounts to the replacement Lt —>■ Lt -I- 2 si = Lt — 2 consistently 
in all previous expressions (6.15), (6.16), and (6.17), as could be easily traced by consistently introducing —m'^ si into 
Eqs. (6.10), (6.12). Therefore, it means that at one-loop RGOPT order si can be simply absorbed by a change of 
scale (or renormalization scheme) definition, /i —>■ Thus apart from changing the reference scale with 

respect to the MS-scheme, it does not really change physical results: if redefining accordingly the coupling with a RG 
evolution /r —>■ we obtain strictly identical results. With this in mind, for now on we proceed with the simplest 

choice Si = 0 at one-loop. 

Goming back to Eq.(6.16) it is obviously an expansion solely in the single parameter ( 6 oA). Using the bo value from 
Eq. (3.4) we have: TTy'2bo/3 = l/(2-\/2), which tells that the first order coefficient differs from the standard Debye 
mass in (2.8), m|) ~ (A/24)T^, being -Ji larger. This originates directly from the correct value bo = 3 /( 167 r^) used in 
a = lo/bo in (4.1), which is the only value compatible with RG invariance. The factor 3 in bo is the statistical factor 
originating from three similar graphs contributing to the /3 function, as is well-known. Thus, the first perturbative 
coefficient of the Debye screening mass would be obtained from Eq. (6.16) if one would take bo/3 = l/( 167 r^) as given 
by a single loop contributing to the self-energy at one-loop, as argued in Ref. [21]. The standard perturbative term 
of order A, comparing with Eq. (2.8), is also reproduced provided again that one takes bo = l/( 167 r^) in Eq. (6.16). 
But, as mentioned above, Eq. (6.16)-(6.17) reproduce exactly (at arbitrary orders) the large Wresults (e.g. Eq. (5.7) 
of ref. [ 22 ]), as can be checked upon identifying the correct large-value of bo = l/( 167 r^) in the normalization of 
[22]. This factor 3 discrepancy in the optimized mass for the N = 1 (j>^ model from bo mismatch is not a problem, 
since the OPT nonperturbative variational mass, not being a physical parameter, has no physical connection with the 
perturbative physical screening mass in Eq. (2.8) and is therefore not required to reproduce the latter. Indeed, there 
are ln(A) terms appearing at the three-loop order in the genuine screening mass[3, 7], Eq. (2.8), that are not present 
in the expansion of Eq. (6.16) which only involves A and A^/^ powers. Incidentally the fact that the standard OPT 
or SPT/HTLpt correctly reproduces the hrst two orders of the thermal perturbative mass expansion [ 6 , 8 ] appears in 
retrospect merely accidental, due to the common canonical choice a = 1/2 in (4.1), i.e. as if one had taken a = 70/60 
with bo —t bo/Sf 

The (exact) expression given by Eq. (6.16) is plotted as a function of the coupling in Fig. 1 in the very common 
normalization[3] A = 24(/^, where it is compared to the standard purely perturbative thermal screening mass, tod/T, 
with scale dependence illustrations. In particular, we remark the saturation of the optimized mass for sufficiently 
large coupling, which agrees qualitatively well with what is obtained at the two-loop order in [21]. This saturation 
can be seen more explicitly by expanding Eq. (6.16) for strong coupling: 




1 -l- |Tt ~ 1 

Lj' 


o{x- 


(6.18) 


The above relation reveals that even if we do not expect our approximation to be valid for arbitrarily large coupling 
the relation fh/T < 1 is always valid while m/T <C 1 in the more perturbative range (see Fig.l). Therefore, the 
high- T approximation used to derive those analytic expressions is fully justified a posteriori. Goncerning the scale 
dependence, in order to compare with standard results, we use the physical OPT solution, Eq. (6.16), replacing 
A = A(^) by its “exact” one-loop running coupling in Eq. (3.9). We take fxo = 27rT as a reference scale, and vary as 
usual the scale /i in the range [ttT, 47rT]. The plot in Fig. 1 is only made for illustration and comparison with the 
standard perturbative thermal mass, since as explained earlier Eq. (6.16) is exactly scale-invariant, i.e. the RGOPT 
mass gap in Fig. 1 has actually zero thickness, being valid for any scale. ® Even if one would use the approximate 
one-loop expanded running coupling (i.e. the last expression in the RHS of Eq. (3.9)), the scale dependence would be 
extremely moderate, barely visible on the same plot. 


We anticipate that for QCD, it happens accidentally that 7 ^(QCD) = bo(QCD)/2, where 7 ^(QCD) is the gluon anomalous mass 
dimension easily calculable from the relevant counterterm given e.g. in [ 8 , 9]. Thus a = 1/2 for the gluonic contributions so that the 
analoguous RG-compatible OPT mass (6.17) will coincide for the first few perturbative terms with the QCD gluon screening mass [3]. 

® Except for the fact that at some large scale p, depending on value, one hits the naive (one-loop) Landau pole, more precisely at 

t/{2ttT) = e^dto^CZirT)). £qj, instance for g = t/A(27rT)/24 = 1 the Landau pole is reached at Tl{2'nT) ~ 8.96. 
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FIG. 1. RGOPT mass m/T at one-loop (<5°) order (thick line), versus standard one-loop perturbative mass as functions of 
gr(/ro = 27rT) = (A(^o)/24)^/^ with scale-dependence. Grey (light blue) bands: variation of the standard perturbative mass 
between jj, = ttT and inT using the exact one-loop running coupling in Eq. (3.9). NB: the RGOPT mass has actually zero 
thickness since it is exactly scale-invariant. 


Next, coming to the pressure it has a rather simple expression at this one-loop order, in terms of the OPT mass rh 
and normalized to the ideal gas pressure Pq = 7r^r"^/90 (still keeping the general 6o-dependence): 


^ 15 


15 rfi^ 45/1 

2t^3 j'S d" \6oA ^ 





(6.19) 


where we remind that this is actually an approximation according to using (2.7), which we argue is however precise 
at the 10“^ level up to x = m/T < 1. One can thus plug the OPT mass expression, Eq. (6.16), into Eq. (6.19) to 
obtain the full A-dependence. After some algebra it takes a compact form: 




Gll + ^G 


3/2 


0{x^) , 


( 6 . 20 ) 


where we defined 1/G = 1/(6oA(t)) +At = 1/(6oA(to))+ 7 B — In2, to emphasize that P is exactly scale-invariant and 
only depends on the single parameter X{fio). Expression (6.20) also explicitly separates the “perturbative” first three 
terms from the clearly “non-perturbative” last term, and is valid implicitly in the high-T approximation as indicated, 
but very precise as long as x < 1, corresponding to G < 3.3 meaning very strong coupling for A(/io) ^ G/bo. The 
first neglected term in (6.20) is actually —15C[3]/(1287r®) x® ~ —1.4710“"^ x® (which is indeed exactly the last term of 
Eq. (5.8) in [22] in the large-iV normalization of bo). One should be evidently cautious not to use (6.20) beyond its 
range of validity, typically for m T, in such a case one rather solves the exact mass gap, Eq. (6.12)). 

One may also easily derive the perturbative expansion of the pressure, which reads to few first orders 

^ ~ 1 - fa + + |(Lt - 6)a^ - §^6 (Lt - |) - | [Lt {Lt - 12) + 6] 

-b [20Lr (At - 3) -b 9] -b |Lt [At {Lt - 18) + 18] + O(a^G) ^ (6.21) 

where a = boX. We remark again that the one- and two-loop standard perturbative terms for the physical massless 
pressure [3, 7, 44] would be reproduced by Eq. (6.21) if effectively replacing, following ref. [21], bo bo/3 = l/(167r^). 
But it is easily seen that the correct coefficients of the leading logarithm terms of the pressure at arbitrary order 
n: a"/n"“^(/t) and lii”“^(/x) (appearing first at three-loop order A^), are given by (6oA/3)(6o ln(/i) A)”“^ 

with the correct bo (comparing with Ref. [44] where the calculation was performed up to order In a using RG 
techniques). Thus taking bo/3 is not consistent beyond the first two perturbative terms. Eq. (6.21) being RG 
invariant by construction correctly reproduces the leading logarithm structure to all orders beyond two-loop order. 
We stress also that Eq. (6.21) correctly reproduces all perturbative terms of the large-result in [22], when taking 
the correct value of the large-TV bo- Since our expressions (6.19), (6.20) are valid for arbitrary N we can in principle 
follow continuously the pressure from large N to N = 1, and while doing this there is no reason to abruptly modify 
the correct bo{N) to some other “effective” bo value. It is useful at this stage to compare this behavior with the SPT 
pressure up to two-loop or higher order [7], basically build on taking a = 1/2 in (4.1): it does reproduce the coefficients 
of the standard perturbative pressure up to second order, but not the correct leading logarithm coefficients at 
order and beyond, as a consequence of missing RG invariance (it would need to rescale In /i —3 In /r to reproduce 










16 


those logarithms). So, the scale dependence of the SPT pressure is unmatched at order c? and beyond if using the 
standard running coupling with bo- 

Thus, keeping the correct bo, which is compelling in our RG-based approach, the optimized results Eqs. (6.20) and 
( 6 . 21 ) differ from the first two terms of the standard (massless) perturbative pressure for small coupling values by 
A(/j,o) —>■ Apert(Mo)/3, for the canonically normalized coupling of the N = 1 scalar model. But this is not a problem, 
simply a different calibration: the exactly scale-invariant RGOPT pressure (6.20) only depends on the single coupling 
G or equivalently X{p,o), still an arbitrary parameter at this stage, since the model is not fully specified by any data 
fixing a physical input scale, /tq. So the pressure as a function of this coupling has a limited physical meaning. 
When going to higher loops, since the nonperturbative RGOPT approximations resum more higher orders, it is not 
surprising that they differ from standard perturbation when expressed in terms of the original perturbative coupling 
A(/r), and one expects to obtain a better approximation for large coupling. The only mandatory feature of any such 
approximation is certainly the Stefan-Boltzmann limit P —>■ Pq for A —0, trivially fulfilled by (6.20). Indeed, those 
features do not contradict truly physical results, as this apparent discrepancy disappears if expressing the pressure in 
terms of the physical mass: to see it, we solve Eq. (4.2) now reciprocally, for A(m), and replacing it in (6.19). It gives 
simply: 


p(i)/Po = 1 - 15a:V(87r2) -t 15a:V(87r3) -p 0(10-V). 


( 6 . 22 ) 


But here x = m/T is arbitrary as we already used (6.15) to fix A(m). Now, taking for m the physical screening 
mass [27] tod in Eq. (2.8), as easily checked it exactly reproduces the first two terms of the standard physical 
pressure [3]. Thus if we would plot our results in terms of the screening mass, P(m|)/T^), we would have very good 
agreement with standard results for sufficiently small screening mass and deviations for larger m?, but the study 
of scale dependence which is our main concern, specially for large coupling values, would be much more difficult. In 
the sequel we keep the results (6.20) in terms of the MS coupling A(/io), which scale dependence is well-defined at a 
given perturbative order, since our aim is mainly to compare the scale dependence with other results in the literature 
also mostly expressed in terms of the running coupling. 

The exact expression for P'XX/Po, Eq. (6.19), is plotted in Fig. 2 where the RGOPT result is compared with the 
standard perturbative expansion at order A ^ g^. The pressure for the rescaled coupling A —A/3 is also shown on 
the same Figure just for the sake of illustration, which accordingly compares better with the standard perturbative 
pressure for small coupling values. The improvement of scale (in)dependence of RGOPT is once again drastic at this 
one-loop order: using the exact one-loop resummed coupling, Eq. (3.9), the RGOPT pressure is exactly scale-invariant, 
which is obvious in Eq. (6.19) since fh is itself exactly scale-invariant and the combination l/A(/r) -|- boLx too, as 
discussed above. This feature is well illustrated by Fig. 2 by comparing the RGOPT with the standard perturbative 
pressure at one-loop which has a notoriously large scale dependence. 

To conclude this section we stress that all the previous RGOPT one-loop results, reproducing among other things 
the exact large-N results, only rely so far on the very simple massive one-loop free energy graph and the knowledge 
of the first order RG coefficients bo, 79 . But these results are not too surprising since the RG properties, if fully 
exploited, involve informations on all orders “daisy” and “super-daisy” foam graphs like those explicitly resummed in 
the large-limit in [ 22 ]. 



rescaled RGOPT 1-loop 
RGOPT 1-loop 
pert. 1-loop 


FIG. 2. One-loop (5°) RGOPT pressure (thick), and rescaling A ^ A/3 (dashed), versus standard one-loop (dotted, light blue 
bands) pressure as function of <? = -\/A(^)/24 with scale-dependence between g = nT and g — AnT. NB: the RGOPT pressure 
has actually zero thickness since it is exactly scale-invariant. 
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VII. 2-LOOP, 0{S) 


We now switch to the two-loop order, thus incorporating all the terms in Eq. (2.4), adding the subtraction terms in 
Eqs. (3.7) and (3.8), and then performing the 5-expansion consistently to order 5 before setting 5=1. Eor simplicity 
we first consider the minimal prescription taking only si 0 , and will also later consider the relevant changes if the $2 
term, formally of order A, is also included. The main novelty at two-loop order is that now the RG relation, Eq. (5.2), 
gives in general a nontrivial additional constraint, that can be used alternatively to the OPT equation, or combined 
with the latter to completely fix m and A in terms of the only remaining free parameter, /i/T (apart from the overall 
dimensional dependence in T). This is another difference with standard OPT or SPT/HTLpt, in which the coupling 
remains undetermined and the generally adopted prescription is to take its perturbative value as a function of /x/T 
and a reference coupling A(/xo) value. In our case, as explained previously one of the RG solutions is matching this 
standard perturbative behavior for A —>■ 0, but for moderate or larger coupling values it will give a nonperturbative 
dependence. One may thus follow at this stage different possible prescriptions: either, one can use any of the two 
OPT and RG equations, to be solved customarily for m(A), next using a two-loop order running coupling, in order to 
compare with other resummation methods. Alternatively, one may consider the full RG and OPT combined solutions. 
Whatever way, neither the optimized mass m solution of Eq. (4.2), nor the optimized coupling A when combining 
the former with the constraint given by Eq. (5.2), have intrinsic universal physical meaning. Both should better be 
viewed as intermediate stage values, to be used only within the physical quantities such as the pressure P{rh, A, ^,/T). 
In particular there is no contradiction between the “fixed” optimized coupling A and the standard running coupling, 
obtained from a different (standard) perturbative RG equation. 


Like for the one-loop approximation, we can express all two-loop RGOPT results in terms of the one-loop self-energy 
defined in Eq. ( 6 .II). After some algebra, the 0(5), 5 —>■ 1 free energy takes a compact form: 


J^o = - 


m 

(4^ 


1 Si . 

A 

3ooA 3 


'p,R 


ln(p2 


/ 270 A 

A V ^0 y 




2A 


2 

R ^ 


(7.1) 


where the index “i?” in the integration means taking the finite part of this already renormalized expression. We also 
kept as much as possible a general dependence on RG coefficients. The first three terms originate from the subtraction 
terms Si in Eqs. (3.7)-(3.8). Notice also the different coefficient l/(36o) as a result of expanding to 0(5^), instead 
of previous l/(25o) at one-loop 5*^ order. As already mentioned we first consider for simplicity S 2 = 0 in the sequel, 
while the effects from S 2 0 (that incorporates a RG part of the three-loop contributions) will be discussed later. 
Next after straightforward manipulations the OPT Eq. (4.2), and the reduced RG operator, Eq. (5.2), acquires a 
compact neat form: 


foprim, A 




+ + So 



= 0 , 


(7.2) 


/rg(w,A, ^) = h 




+ i/3(2)(A)A2 = 0, 


(7.3) 


with h = ( 47 r)“^, + biX^ is the standard /3-function restricted to two-loops, and recalling also that 

Si = —1. We have defined for convenience the reduced (dimensionless) self-energy S{m,fi,T) = Efl/(m^A) thus 
independent of A, which makes the coupling dependence very transparent in Eqs. (7.2) and (7.3). For completeness 
and further use below we also have from Eq.( 6 .II): 


Note that in principle the reduced RG Eq.(7.3) is only valid when combined with Eq. (7.2) since the latter removes 
the dm part of the complete RG operator in Eq. (5.1). Glearly, in the above normalization, the complete RG Eq. (5.1) 
reads 


/RGfull = /rG + 27m(A)/oPT = 0 , (7.5) 

where the anomalous mass dimension 7 ^ (truncated at the two-loop order) was defined in Eq. (3.3). Therefore, to 
obtain the most general solution to(A) consistent with arbitrary coupling values, Eq. (7.5) should be solved. As we 
shall see below, the solutions torg(A) and mR,Gfuii(A) are very close for sufficiently small A but can depart substantially 
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from each other for arbitrarily large coupling. The advantage of using the reduced RG operator is that the solution 
can be more easily found when looking for the intersection between the two RG and OPT solutions toopt('^) and 
torg(A). 

Before proceeding, another digressing remark is that the standard OPT/SPT would correspond to a much simpler 
OPT equation than Eq. (7.2), since in particular the first two (subtraction) terms soiSi would be absent, resulting 
in a simple OPT self-consistent solution: = Y,r. Moreover the modified RG Eqs. (7.3) or (7.5) are usually never 

considered within the standard OPT/SPT or HTLpt applications. The fact that the OPT and RG relations, Eq. (7.2) 
and Eq. (7.3), are more involved is expected since the one-loop RGOPT already gives nontrivial results qualitatively 
similar to two-loop standard OPT/SPT. Accordingly, the relative complexity of RGOPT equations at two-loop is 
due to the more information they carry on higher RG orders, and is a price for a more efficient and RG-consistent 
resummation procedure. At this stage the OPT and RG Eqs. (7.2), (7.3) could be solved exactly for X{m, being 
respectively quadratic and cubic algebraic equations in A. But to compare with most results in the literature it is more 
customary to rather solve for a mass gap m(A, /i/T), to obtain in a next stage the pressure or other thermodynamical 
quantities as a function of the coupling. 

Considering first the (reduced) RG Eq. (7.3), it is best solved in a first stage as a simple quadratic equation for 
S'(m,/i, T) whose mass gap solutions are: 


'^rg(A) — 


Sr 

Am^ 


iTyr^aTfAKTr^ 

3A(l + tA) 


(7.6) 


where again the explicit RG dependence has been kept for generality. Note that the coupling dependence is entirely 
contained in the RHS of Eq. (7.6) since by definition S does not depend on A, prior to using the RG Eq. (7.3). Just to 
see where we stand before considering the more involved exact mass gap solutions at two-loop order, let us consider 
Eq. (7.6) by crudely neglecting the two-loop /3-function coefficient, bi = 0, or equivalently taking the leading term when 
expanding (7.6) for A —>■ 0 . It gives immediately S'rg(A) = 1/(3 A) which, recalling that S{m,n,T) = ^^/(Am^) 
and using (6.11), is nothing but the one-loop mass gap Eq. ( 6 . 8 ) consistently recovered. Once we use a nonzero 
6 i ^ 0 in Eq. (7.6), at it should at two-loop RG order, one could consider a perturbative expansion of Eq. (7.6) to 
gradually include higher perturbative order corrections to the one-loop solution ( 6 . 8 ), but it is algebraically not more 
complicated to solve Eq. (7.6) exactly. Now in fact Eq. (7.6) also reflects a possible complication appearing at two-loop 
order for large coupling values, due to the two-loop ultraviolet fixed point (UVFP) at A = —bo/bi since bi < 0, see 
Eq. (3.4). This purely perturbative UVFP is totally spurious, not only since it disappears at three-loop level (where 
the next coefficient 62 is positive[28] and large enough so that possible non-trivial fixed points are complex) but more 
generally since the existence of nonperturbative UVFP is excluded by the numerical evidence for the triviality [47] 
of the theory. Nevertheless, since the RGOPT construction basically relies on perturbative RG properties, one 
may worry that some of our two-loop results could be affected, if driven by this spurious UVFP. Indeed the (-b) 
solution in Eq. (7.6) is singular at the fixed point value of A, which means that to^/Sr(to^) —>■ 0 i.e. m —>■ 0, while 
the (—) solution is regular: *S'[/g(—& 0 /& 1 ) ~ &i/( 6 &o), which means that it is a priori the solution not wrongly driven 
by the UVFP. Thus, in the sequel we should be careful to identify any behavior that could be an artifact of this 
perturbative fixed point. Now, in terms of the rescaled coupling A = 2Ag'^, the UVFP is at ^ 1.866, and the 
maximum of the /3-function (beyond which the coupling is driven to a really wrong behavior), is at A = —26o/(3&i) 
i.e. g ~ 1.524. Both values are to be considered very large couplings, where the validity of a resummation procedure 
is anyhow questionable. Therefore, as long as one stays safely below these large coupling values, say not too much 
above ~ 1 in practice, our results should remain valid. Moreover, here we are basically focusing on the RG/scale 
invariance issues in very general terms, rather than on the peculiar nonperturbative dynamics of the 4>^ model, which 
is beyond the present scope. Another property to notice is that the exact RG solutions, Eq. (7.6), become complex 
at a coupling Ac = —36o/(2&i) = (47r)^27/34 corresponding to ~ 2.28, thus irrelevant since located beyond the 
fixed point anyway. Yet one should keep in mind that independently of the presence of non-trivial perturbative fixed 
points, complex optimized RG solutions are unavoidably expected to occur at some higher perturbative order from 
exactly solving the OPT and RG equations, as discussed above. 


A. T = 0 

To get some more feeling we first explore the T = 0 case, which is much simpler sin ce bo th RG and OPT equations 
can be solved analytically, e.g. for lnm//x in terms of the (rescaled) coupling g = ^ A/24 (or its reciprocal function 
A(lnTO/^)). The (real parts) of the OPT and RG exact solutions, expressed as h\[m{g )/are plotted in Fig. 3. We 
also show for comparison the one-loop solution (dotted) from Eq. (6.7). The physical branch solutions are clearly 
identified, i.e. those matching standard perturbation for /x <C to with fixed to: Inm/jj, —>■ -boo for A —>■ 0+ (see 
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Re(Ln(m/ (J,)) 



FIG. 3. T = 0 two-loop (5^) order OPT (dashed) and RG (thick) (real parts of) solutions ln^{g), g = ^yX/24:). Also shown 
is the one-loop solution (dotted) from Eq. (6.7). 


also the discussion after Eq. ( 6 . 8 )). One also sees the asymptotic value \nm/fi ~ 1/2 reached for large coupling 
g consistently with Eq. (6.7). The two-loop OPT (dashed) physical branch becomes complex for g > 2.094 (where 
there is a corresponding bifurcation on the figure), in fact very close above the intersection between the two (real) 
physical branch solutions, occurring at 5 ~ 2.08,ln(TO/^) ~ 0.083, thus at fh/g close to 1 and a quite strong coupling 
value. The OPT branch is real again at about 5 > 6 (where the dashed curve shows a little bump). The RG physical 
branch becomes also complex at a slightly higher g ~ 2.28 value, as already noted above after Eq.(7.6). There are 
also two other combined RG and OPT solutions (intersections) sitting on the complex branches, see Pig. 3, but these 
are to be considered unphysical solutions since not connected with the perturbative branches. One can also note 
that the RG br anch has a p ole behavior around 5 — 2, which is a consequence of the above discussed perturbative 
UVFP at 5 = -\/—6o/(246i) ~ 1.866.® The reciprocal function giXnm/g) would show a frozen behavior at 5 ~ 2. In 
fact, we stress that the behavior around 5 ~ 2 , including the solutions becoming complex, is all driven by this naive 
perturbative two-loop UVFP, so that one should simply not trust what happens for g close to those values, say for 
g > 1 to be on the safe side. (Notice however that on the figure the one-loop-like behavior is recovered for much 
larger g values). In particular, the above mentioned real OPT and RG intersection solution at 5 ~ 2.08 is beyond 
the UVFP, thus very untrustable. From examples in other theories [17], we expect that at higher orders the RGOPT 
intersection solution may decrease below the UVFP and stabilize to a more reasonably perturbative value. 

Switching on the thermal contributions modifies coefficients of the relevant RG and OPT equations, which will result 
in real intersecting points for the RG and OPT solutions with somewhat lower coupling values 5^1 for generic g/T 
values as we examine in next subsection below^. 


B. T / 0 


Considering now the thermal contributions, one may solve numerically Eq. (7.6) (or the full RG Eq. (7.5)) and the 
OPT Eq.(7.2), using the exact expression S{m,T,fi) from Eq. (6.11), to obtain x = rfijT as function of \{g) at some 
chosen scale g. Concretely, to solve the OPT gap-equation exactly for arbitrary temperature, it is convenient to first 
solve Eq. (7.2) as a linear equation for S{m/T, fi/T) in terms of A and S'{m/T, g/T), giving trivially 


S = 




1 


(7.7) 


87r2 (2 -b 3A S') ’ 

to be then solved numerically as a mass gap rh{X,g/T) as function of the coupling and scale, using the expressions 
of S, S' in Eqs (6.11), (7.4). In the right hand side of Eq. (7.7), taking only the first term, dominant for A —>■ 0, one 
immediately recovers again the one-loop RGOPT mass gap solution Eq. ( 6 . 8 ) (just like for the above discussed RG 
mass gap (7.6) when taking 61 = 0), while the second term of Eq. (7.7) clearly gives higher order corrections if seen as 
a perturbative expansion. (However, for large coupling values of order g = y^A/24 ~ 1 such a perturbative expansion 
of Eq. (7.7) would not give a very accurate mass gap solution so it is better to solve it exactly numerically). 

The exact OPT Eq. (7.7) and RG Eq. (7.6) are illustrated with their roots for a rather strong coupling value 5 = 1 


® The pole of \nm/g is not exactly at the fixed point of /3(A) due to the mass-dependence entering the RG Eq. (5.2). 

^ When complex solutions occur on physical branches, one may recover real solutions by performing a perturbative scheme change, as 
done in [17]. But this more involved course of action can be avoided in the 0“^ case, at least at the two-loop level. We anticipate however 
that for thermal QCD, RGOPT will unavoidably gives complex solutions, mainly due to the opposite signs of the bo, 61 coefficients due 
to asymptotic freedom. 
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in Fig. 4. As one can see, in general both the two-loop order OPT and RG equations have three real solutions, until 
two solutions become complex (conjugates), which happens for g > 2.09 and g > 2.28 respectively for the OPT and 
RG equations. Since this is well beyond the fake perturbative UVFP, we cannot trust the detailed consequences near 
such large coupling values. For more moderate coupling values as illustrated in Fig. 4, there is one RG and one OPT 
solution with very large m/T ^ 1, in fact behaving for small A as to ^ gi/l^'oA) g^g gg^ easily trace even from the 
one-loop mass gap Eq. (6.12): for to T the T^Ji{m/T) term in (6.12) becomes negligible and one simply recovers 
the T = 0 solution for to(A). But for T 0 this solution does not have the property of a thermal mass, to —>■ 0 
for T —>■ 0, so that the other roots are the physically relevant ones. The behavior of the two intermediate and large 



FIG. 4. Roots of the two-loop (d^) exact OPT Eq. (7.2) and RG Eq. (7.3), as compared with one-loop OPT Eq. (6.12) for 
g = V^(A‘o)/24 = 1, go = 2tvT. The y-axis values are those of /opt Eq. (7.2), /rg Eq. (7.3), and Eq. (6.12) in convenient 
common units. 

values roots is qualitatively similar to the one-loop OPT also illustrated and large-mass gap solution [22], indeed 
recovered as above explained at one-loop order. What is new as compared to one-loop order and seen on Fig. 4 are 
the two extra roots with the lowest x = m/T values both for the OPT and RG equations. Goncerning the RG root 
with the lowest x ~ 0.7 value, it is easily traced to be the one driven by the UV fixed point: for A — —bo/bi, it gives 
a: —?> 0, and we should reject it accordingly. The other OPT root with lowest x ~ 0.6 value is more special: it is 
more easily analyzed by solving the exact OPT Eq. (7.2) for A(to/T), which is a simple quadratic equation. It can 
also be seen to correspond to the perturbatively odd situation where the second term in Eq. (7.7) dominates over the 
first term. Then matching with the perturbative behavior, it gives an “ultrasoft” mass m/T ~ (7r6o/2)A + 0{\^) for 
A —>■ 0, which contradicts the expected behavior of a thermal mass to ^ V^T on general grounds, that was indeed 
found at one-loop order. At higher orders this optimized mass has a standard power series in A. So we consider this 
solution as a spurious unphysical one, an artifact of the more involved two-loop OPT equation. Therefore at two-loop 
order we identify unique physical OPT and RG solution, which are real and positive for all relevant coupling values 
and compatible with the perturbative behavior of a thermal mass, corresponding to the two intermediate value roots 
near a: ~ 1 illustrated on Fig. 4. Once correctly identified, those physical solutions have a qualitative behavior not 
drastically different from the one-loop order OPT solution of (6.12), apart from the more involved algebra. Note also 
that solving Eqs. (7.6) or (7.7) for A(to, ln(/i/r)) one can check analytically that for fixed m/T values, the physical 
A(/i/T) solution decreases logarithmically in the two extreme limits g <C 2ttT and g ^ 2'kT, having a maximum in 
between, a behavior qualitatively quite similar to the large-Af or one-loop case. 

Alternatively for any solutions, as long as a; < 1, one can solve analytically both RG and OPT equations but using 
the high-T approximation (2.7) and derivatives for the relevant thermal integrals, giving 

S(xsm/r<l)=.-^(tT + ^-^) . (7.8) 

In this way one obtains respectively quartic RG and cubic OPT equations in a: = m/T, which are not particularly 
telling but gives algebraic solutions. We checked that all the approximate high-T RG solutions for to and the pressure 
P/Pq for relevant scale values ttT < g < AttT are excellent, departing below 0.1% from the exact T-dependent 
solutions, at least up to a value of the (rescaled) coupling g = y^A/24 ^ 1.5, simply because x remains always lower 
than about ^ 1. Goncerning the OPT solutions, they can give a; > 1 at large coupling g > 1, particularly for the 
higher scale choice g = 47rT, in which case the high-T approximation starts to fail and we better use the exact 
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T-dependent numerical solutions ® of Eq. (7.7). 


1. Comparison with standard perturbation theory 


The OPT physical solution has the following perturbative expansion 

^ ^ - 2Tt)A3/^ - ^^(3 - 4Lt)A2 + • • • (7.9) 

consistently with the first two terms of Eq. (6.17). Concerning the terms of order A^/^ and A^, the (leading) logarithms 
coefficient of Lp are the same with respect to the one-loop expansion Eq. (6.17), as it should be to all orders. The 
differences appear only in the constant terms, due both to si 0 and other terms of order A*^ in the original expression 
of the free energy (7.1). Note that the relative weight of Si is changed at order 5^, with si/3 in (7.1). Quite similarly, 
incorporating S 2 7 ^ 0 gives differences only visible at order A^ not given in Eq. (7.9), so it has very little effect at least 
on the perturbative mass expansion. Plugging this in the pressure expression and expanding perturbatively (while 
keeping the coupling free for the moment) we obtain the ratio of the two-loop RGOPT pressure to the ideal gas one 
Po = 7r2r4/90 as: 


P, 


( 2 ) 


OPT 

Pn 


l-^a-b^V6a^/^-b^(PT-20/3)a^-^y6(Lr-13/6)a®/2-^[LT(PT-40/3)-f44/3]a^-fC>(A^/2) , (7.10) 


where a = 6oA, giving the same first two perturbative terms as the RGOPT one-loop result Eq. (6.21). We remark 
also, just like in the one-loop case, that if rather solving either the OPT or the RG Eqs. (7.2), (7.5) reciprocally for 
A(m), and then taking for m the physical screening mass Eq. (2.8), one consistently recovers the standard perturbative 
expression of the physical pressure up to two-loop order. Remark that the effect of incorporating S 2 7 ^ 0 instead of 
the simplest choice with only si, appears first in the very last non-logarithmic term of order m'^A ^ A^ in (7.10), 
replacing 44/3—7 637/45 — 4C[3]/5, which has very little effects upon the final numerical results as long as A remains 
reasonably perturbative. This is completely expected since perturbatively m‘^S 2 X in (3.7) is of lowest perturbative 
order ^ A^S 2 A ^ S 2 A^. 

Alternatively, we also consider the solution m{X,fi/T) obtained from the (full) RG Eqs. (7.3) and (7.5). The latter 
equations can be solved for the exact T-dependence first as quadratic equation for S{m/T, fi/T) and then numerically 
for m{X, fi/T). But as mentioned above the high-T approximation is excellent in this case for all relevant coupling 
values. The high-T approximated RG equations gives two negative and two positive rh/T solutions, the latter having 
both the correct perturbative behavior, but with one saturating faster for large coupling, which is illustrated as 
“solution 2” in Fig. 5. In Fig. 5 we also show the scale dependence in the range ttT < /r < 47rT for all solutions 


m/T 



FIG. 5. Two-loop (J^) optimized mass solutions obtained from OPT Eq.(4.2) (dashed, light blue on line) and RG Eq.(5.2) 
two solutions: thick lines (green on line) and dot-dashed (red on-line) respectively, as function of g = ■y/A(/r)/24 with scale 
dependence ttT < p < 4nT. 


® Such remarks are important to keep in mind in view of possible applications to thermal QCD, for which beyond the one-loop level only 
high-T expansion results are available [9]. 
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FIG. 6. RGOPT one-loop (black line) and two-loop (dot-dashed, green on line) versus standard perturbative pressure P/Po{g = 
\/A/24) at one-loop and two-loop. The different bands give the scale dependence between fj, — ttT and fi = dvrT. NB: the 
RGOPT one-loop curve has actually zero thickness since it is exactly scale-invariant. 


X = rh{g, ii)/T. The “lower” values RG solution 2 happens to reach a maximum and then decreases towards zero 
for large coupling g ^ 1.86 (not shown in the figure), an odd behavior which is in fact completely driven by the 
perturbative two-loop UVFP. Incidentally this solution exhibits an extremely small, almost totally negligible scale 
dependence up to 5 ^ 1, as can be seen in Fig. 5, in agreement with what is intuitively expected near a RG fixed point. 
The corresponding pressure, if plugging this solution within (7.1), has even smaller scale dependence. Although we 
shall simply discard this solution, since the naive two-loop UVFP contradicts the genuine nonperturbative dynamics 
of the (j)^ model, it is worth remarking that the corresponding RGOPT solution faithfully reflects the quasi scale 
invariant behavior of the UVFP. The remaining (then unique physical) RG solution 1 exhibits a more pronounced 
scale dependence as seen in Fig. 5, which will be further explained below. It has a perturbative expansion for small 
coupling with the first two order terms identical with the other OPT solution above, Eq.(7.9): 

# - A3/^ + 5^(5 + lOLriA^ + ... (7.11, 

where one can easily check that the coefficients of Lt, thus of In/r, are identical, which should be the case to all orders 
consistently with RG invariance properties. 

The perturbative expansion of P/Pq is identical to (7.10) for the first few orders and for Lt coefficients at all orders 
if we use instead the rh{X)/T solution, Eq. (7.11). 

Remark in Eig. 5 the intersections, for different fi scale values, between the RG and OPT solutions, which only exist 
for the higher of the two RG solution branches (namely the physical solution not influenced by the perturbative UV 
fixed point). Those intersections can be considered as the full RGOPT solution, which is unique for a given fj, input 
scale. We shall come back later on those full RGOPT solutions, while for the sake of comparison with standard 
perturbation results we consider at the moment the RG and OPT solutions separately, as given functions of the 

coupling. _ 

In Fig. 6 we plot the exact two-loop RGOPT pressure P/Po{g = y/X/24) (dashed lines), as obtained from the 
mass solution of the full RG Eq. (7.5), compared with one-loop RGOPT and standard perturbative one- and two-loop 
results, with scale dependence in the range ttT < /i < 47 rT. To study the scale dependence we use a standard running 
coupling exact at two-loop orders. The generic exact expression for the two-loop running A(/i), generalizing the exact 
one-loop running in Eq. (3.9), can be expressed (see e.g. Ref. [31]) in terms of the (implicit) Lambert “function” [45] 
W{x) = ln(a;/lU): 
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where L = ln/i//ro and where we consider as usual the reference scale /io = 2'kT. Actually Eq. (7.12) gives no visible 
difference (at least up to the relevant moderately large coupling values here studied and the moderate range of scale 
variation) with a more handy perturbatively truncated expansion at order A^: 

hlhiL^^ + 0{X^) . (7.14) 


A ^(/x) ~ A — boL — {biL)X — —b^biL^ 


2t2 


y - 1 ^biL 


The RGOPT improvement on scale dependence with respect to standard perturbative results is drastic for the pressure, 
as one can clearly see in figure 6, although scale invariance is not exact at two-loop like it is at one-loop order: there is 
a moderate residual scale dependence, clearly visible in Fig. 6 for moderately large (rescaled) coupling values g > 0.6. 
What is also clearly seen in Fig. 6 is the much better stability of the RGOPT results, since up to ~ .5, both the 
scale dependence and the difference between one- and two-loop RGOPT pressure are hardly visible at the figure scale, 
in contrast with the already poorly convergent standard perturbative pressure for those values. Accordingly, there 
are important extra cancellations of the scale-dependence happening when m(A, /i) optimized solutions from Fig. 5 
are plugged into the pressure expression P/PQ{m, X, /i) in Eq. (7.1). But we emphasize again that the optimized 
masses mQprp(A,/i) or rh^Q{X, g) at order-5^ are intermediate, unphysical quantities, therefore not expected to be 
themselves scale-invariant in general (although at one-loop order, mQprp(A,/x) in Eq. (6.16) is exactly (one-loop) 
scale-invariant as explained above). Within the RGOPT procedure, RG-invariance is by construction required only 
for the pressure which represents the actual physical observable, resulting in Eqs. (7.3) and (7.2) optimizing the 
pressure. Accordingly, the further cancellation of scale-dependence of the two-loop OPT or RG masses ?7iopt(A, /x) or 
wrg(A, /x), once plugged into the pressure P/Po{m, A, /x), is expected from the resummation properties of the RGOPT 
procedure. 

The residual scale dependence of the two-loop RGOPT pressure is unavoidable due to the RGOPT construction being 
not exact but resulting from the optimization of actually two-loop restricted basic free energy by construction, where 
terms of order A^ and higher are truncated. (In contrast the one-loop results above were exactly scale invariant 
because of the perfect matching of the exact one-loop running coupling with the mass gap (6.16)). 

One can make those statements more precise by studying exactly at which perturbative order the scale dependence 
reappears: examining the perturbative expansion of the pressure above in Eq.(7.10) in which the coupling is replaced 
by its running expression at truncated two-loop order (so only with bo, bi dependence) Eqs. (7.12) and (7.14), we have 
checked explicitly that the leading scale dependence reappears first at order A^: 

= (85 1„^-2190.5)+0(V/»). (7.15) 

i.e. formally at four-loop order, that is one order higher than the naively expected three-loop A^ from standard RG 
invariance properties®. This feature can be anticipated in fact with a little insight without involved explicit high order 
expansion calculations: recall that the RGOPT construction including the subtraction terms in (3.7), together with 
the RG invariance preserving interpolation (4.1), guarantees that the two-loop free energy (2.4) is RG invariant up 
to neglected three-loop terms, of order (!l(m^A®), but for arbitrary mass m. This implies that the mass gap obtained 
either from the RG Eq. (7.6) or OPT Eq. (7.7) has a remnant scale dependence appearing perturbatively at order 

77 j 2 /y 2 ^ A(l-|-1-A® In^), as could easily be checked by explicit expansion too. So it means that the lowest possible 

order at which a remnant scale dependence appears in the free energy (2.4) is given by the first subtraction terms 
—som'^/X, giving remnant scale dependence at perturbative order A^ln/x). In contrast the remnant scale-dependence 
of the standard perturbative two-loop pressure appears at the expected order A®. 

Using alternatively the RG solution m)^Q(A), we find a residual perturbative scale-dependence reappearing also consis¬ 
tently at order A^ with the same coefficient of In /x than in Eq.(7.15), and a very similar constant term, with ^ —2190.5 
in Eq. (7.15) replaced by ~ —2009.9. However, we stress that Eq. (7.15) allows to understand the perturbative be¬ 
havior of the remnant scale dependence, but does not properly reflect the actual nonperturbative scale-dependence of 
the full RG or OPT solutions, which we used for the plots in Fig. 6. Indeed the actual values and scale dependence 
of the pressure for relatively large coupling g > 0.6 are very different than what would be obtained by a finite order 
perturbative truncation at order A^. Accordingly, the moderate residual scale dependence seen on the plots for large 
g values appears much better than what (7.15) would give, as a result of further nonperturbative cancellations among 
successive higher orders (of course higher orders of the RGOPT resummation anyway based on a two-loop truncated 
free energy). Glearly when g becomes of order 1 fh/T is also of order 1 so that the previous perturbative reasoning 


® When rescaling the coupling as A —)■ 2Ag^, the leading remnant scale-dependence in Eq.(7.15) gives: ^ (0.075 In/i/xio — 1.92)g®, where 
the lowest order coefficients are roughly of order 0(1) in this normalization. 
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FIG. 7. Two-loop pressure P/Po{g = iyA/24) obtained from the full RG solution of (7.5) (dot-dashed, green) and obtained 
from the OPT solution (4.2) (dotted, light red), with scale dependence between g — ttT and g — AnT. The RGOPT one-loop 
pressure (black line) is also shown for comparison. 



FIG. 8 . Same as in Fig. 7, but with S 2 7 ^ 0 in Eq. (3.7). 


with (7.15) no longer apply and we simply recover that the scale invariance is guaranteed by construction up to 
remnant terms of order m^}? In fj,. 

As another illustration, in Fig 7 the exact pressure obtained from the OPT mass solution of Eq. (7.2) is compared with 
the one obtained from the RG equation. Although both RG and OPT solutions have perturbative scale-dependence 
reapparing at the same order, like in (7.15), the nonperturbative scale dependence of the OPT solution is almost 
negligible until g ~ 0.8, while for larger coupling values it becomes more important than the one of the RG solution 
(as one can see also for the corresponding OPT mass in Fig. 5). We also remark in Fig 7 that the pressure obtained 
from the RG mass remains closer to the one-loop RGOPT pressure for large coupling values, so that the convergence 
appears better. This, together with the better scale independence for large coupling, is not very surprising since the 
RG mass solution originates from the RG Eq. (7.5) at the two-loop level. In contrast, the OPT solution results solely 
from the mass optimization, Eq. (4.2), which incorporates RG properties more indirectly: it also exhibits good scale 
invariance up to relatively large values g > 0.8, beyond which it degrades quite rapidly. Also, the fact that the RG and 
OPT solutions are very close to each other until relatively large coupling values g ~ 0.6 shows an overall consistency, 
by quantifying the relatively small lack of exact RG invariance, since for an exact nonperturbative result the OPT 
and RG solutions would be identical. 

We now consider incorporating the S 2 7 ^ 0 term from Eq. (3.7), which formally belongs to the two-loop order. 
We already mentioned above that perturbatively it evidently only affects the (re)-expanded pressure (7.10) at order 
rh'^X ^ A^. Accordingly, it does not affect the perturbative order at which the moderate residual scale-dependence first 
reappears, Eq. (7.15). However, the nonperturbative RG and OPT pressures are affected for larger coupling values, 
as intuitively expected since including S 2 7 ^ 0 incorporates a (RG) part of the three-loop contributions. Indeed, as 
shown in Fig. 8 , it improves slightly the (nonperturbative) scale dependence of the RG pressure (long dashed), in 
addition the OPT and RG pressures are closer to each other for large coupling values (but the scale dependence 
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FIG. 9. one- and two-loop RGOPT pressure versus one and two-loop standard perturbative pressure and two-loop SPT pressure, 
with scale-dependence. 


is still better for the RG solution than for the OPT one). More remarkably, with S 2 ^ 0 the two-loop pressure 
obtained from the RG mass gap is seen to almost coincide with the one-loop pressure, up to relatively large coupling 
g ^ 1, thus improving the apparent convergence further more. Actually, this rather spectacular coincidence with 
the one-loop pressure for a large range of coupling values is partly accidental: after applying (4.1), there are some 
partial cancellations of the two-loop contributions happening for relatively large coupling values (due to opposite signs 
si,S 2 ), between si/3 = —1/3 and the dominant two-loop term S 2 A = 2Ag‘^S2 (see (7.1), (3.7), (3.8)), with a maximal 
cancellation for g ~ 1, with 24g^S2 — 1/3. Since the Sk values depend on the particular RG coefficients, in a different 
theory si and S 2 may have the same sign, or rather different magnitudes, possibly giving no such partial cancellations. 
Nevertheless the coincidence with the one-loop pressure is excellent even for relatively large intermediate coupling 
values e.g. g ~ 0.6 — 0.7, where 24g^S2 — .12 — .16 does not make the cancellation much effective. Accordingly, there 
is also clearly a more generic effective stabilization of the perturbative series resulting from the RGOPT construction, 
with an improved convergence and scale dependence when incorporating higher RG order dependence, as intuitively 
expected. 


2. Comparison with other (SPT and HTLpt) variational resummation approaches 

This stability and improved scale-dependence is also illustrated in the next Fig. 9, where the RG pressure is 
compared with a standard two-loop OPT/SPT [7, 14]: more precisely, discarding Sq in (3.7), taking a = 1/2 in (4.1), 
and using solely Eq. (4.2). To compare with another mass prescription, instead of the mass optimization we use the 
screening mass, Eq. (2.8) (consistently truncated at two-loop order), plugged in the expression of the free energy, 
similarly to the prescription mostly adopted for QCD HTLpt [9]. One sees that using the optimized mass within 
the SPT/OPT gives a better scale dependence, although SPT with optimized or screening mass both have a definitely 
stronger scale dependence than the RGOPT for moderately large coupling values. To quantify what is illustrated 
in Fig. 9 more precisely, let us indicate the relative scale variation of the various methods for the relatively large 
(rescaled) coupling value g{2'KT) = 1: the corresponding variation of P/Po between g = ttT and 47rT is ~ 8%, 0.8%, 
1.5%, 0.35%, and 0.3% respectively for the 2-loop standard perturbation, SPT (optimized mass), SPT (screening 
mass), RGOPT, and RGOPT with S 2 ^ 0. Beyond g > 1 the two-loop RGOPT scale-dependence increases more 
rapidly, but is still only ^ 1.610“^ for g ^ 1.5 (while for such large coupling the relative variation of the standard 
perturbative two-loop pressure is as large as ^ 0.25, and the scale variations of the HTLpt and standard OPT methods 
become very important too). It would also be interesting to compare quantitatively our scale variation results with 
the residual scale-dependence appearing in the 2PI approach at three-loop order [23]. However, a precise comparison 
appears difficult, since the renormalisation scheme and scale used in [23] are completely different and not easy to 
translate into the present scale variation in the MS-scheme. 

The RGOPT improvement on scale dependence at the two-loop order may not appear so spectacular as compared 


10 


Using the perturbative screening mass instead of the optimized mass gap is essentially the procedure in HLTpt applications to QCD 
because the optimized mass has no real solutions [9] 
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with SPT, merely a factor ~ 2 — 4 for g ~ 1, essentially because the two-loop order SPT scale dependence here 
illustrated is still moderate. But at the three-loop level the SPT scale dependence becomes much larger [7] (and 
similarly HTLpt QCD[9]), of order 1 for moderately large coupling. As previously explained this can be traced 
both to a ^ 7o/^o in (4.1) together with the missing RG invariance from m^ln/i terms for arbitrary m in (2.4): 
since perturbatively m^ln/i ~ A^ln^, the lack of scale invariance at formally one-loop order in the original massive 
free energy remains somewhat screened at one- and two-loop thermal perturbative expansion order, until it plainly 
resurfaces at three-loop order. Remarkably the SPT pressure has even been calculated more recently to four-loop 
order [46], and these results show very good convergence with respect to two- and thee-loop SPT for the central scale 
choice n = 27rT. But the remnant scale dependence is not illustrated in [46]. According to our general arguments 
we do not see how the missing RG-consistent one-loop terms could be compensated by going to higher orders if not 
present from the beginning. In contrast the RGOPT scale dependence remains very moderate, as illustrated in Fig. 
9 when including a RG part of the three-loop contributions from S 2 0. Moreover without explicitly calculating 
the full three-loop or higher order RGOPT pressure, we can be confident that the scale-dependence should further 
improve at higher orders: being built on perturbative RG invariance at order k for arbitrary m, the RGOPT mass 

gap will exhibit remnant scale dependence as ~ Ar^(l H-1- 0(A^+^ In^)), thus the dominant scale dependence 

in the free energy (namely the lowest perturbative order at which scale dependence will resurface), coming from the 
leading term —som'^/X, should be of order 0{X'^~^'^). 

Goming back to the two-loop order pressure, we have also checked the variation of our results against various pertur¬ 
bative truncations: for instance given that all our calculation relies on the basic two-loop free energy (2.4), it may be 
unnecessary rehnement to use the exact running coupling as in (7.12). Thus we looked for variation when truncating 
the running Eq. (7.14), at order A^, or A^ only. Also, we studied the effect of truncating in the RG Eqs. (7.3) and (7.5) 
terms of order A^, originating from the biX^d\ term in (3.2), as it is formally one order higher than perturbatively 
strictly required (and the 6i-dependence enters anyway at lower orders in the RG Eq. due to the —sq/A subtraction). 
We obtain very good stability, since the maximal resulting variations, for the relatively large coupling g = 1, and 
scale choice g, = 47rT, is below 10“^ for Prg/Pq, and somewhat worst but reasonably below 10“^ for Pqpt/Po- 


3. Full RGOPT two-loop solution 


(o') 

Einally, we can calculate PrgoptZ-^o for the complete two-loop RGOPT solution, given by the non-trivial inter¬ 
section between the RG and OPT equations (as illustrated in Pig. 5), as function of the only free scale parameter, 
that we choose as t = g,/{2nT) (so the standard central scale choice corresponds to t = 1). Ideally for an all order 
calculation with exact scale invariance one would expect both RG and OPT equations to give identical m(A,/i/T) 
solutions. This is indeed the case for simpler models, like typically [15] the large Wlimit of the 0{N) GN-model. But 
for a more involved theory at a finite order one expects some differences between the RG and OPT solutions since RG 
properties are only imposed perturbatively, those remnant differences reflecting the lack of exact RG/scale invariance. 
Just like the stationary mass solution is expected to approximate the actually massless result, the intersection between 
the OPT and RG curves at a given order, defining (variational) “fixed-point” mass and coupling, is expected to give 
a sensible approximation to the exactly scale-invariant nonperturbative result. 

Por the standard central scale choice t = 1, the solution corresponding to the unique intersection of the physical RG 
and OPT branch solutions, readily seen for rh and g in Pig. 5, gives: 


^ 0.907 , 
Pq 


(7.16) 


obtained using the simpler high-T expansion solutions. (NB for t > .5 approximately, x{t) = rh{g,)/T remains smaller 
than 1, justifying a posteriori the high-T approximation. Typically, for x = .9 the high-T approximation is already 
correct at the level of 0.2% and differences are completely negligible for smaller x.) 

The result in (7.16) needs further comments on its physical interpretation. Recall that the truly nonperturbative 
massless pressure expression P{X{g/T), n/T)/P q, if that was available, would actually be a function of the single 
coupling A(/io = SttTq) given at some input scale, thanks to exact scale invariance: for any renormalization scale 
choice fj, the nonperturbative running A(/x) would exactly compensate the explicit /r/T dependence. Incidentally, this 
is precisely the situation happening at the one-loop RGOPT order, where the exactly scale invariant P/Pq in (6.20) 
only depends on the single parameter 5oA(/ro). Now, at higher loop orders in the standard perturbative approaches 
(or similarly in the SPT/HTLpt approaches after expressing the mass gap to(A) in terms of the coupling), due to 
imperfect scale invariance giving remnant perturbative terms ln^(/r/(27rT)), one avoids large logarithms by fixing /x 
of order ^ 27rT, which makes the scale effectively T-dependent and allows to study P/Po{X{T/To)) by varying the 
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coupling as a function of the scale/temperature. In our case, using the second RG constraint Eq. (5.1) enforces 
RG invariance at a limited (here two-loop) perturbative order, mimicking exact RG invariance. At T = 0 two- 
loop order, as already mentioned in Sec. VILA there is also a non-trivial intersecting OPT and RG solution at 
g ~ 2.08, ln(TO//x) ~ 0.08 (see also Fig. 3), that gives the free energy as Jb oc which still requires a reference 
physical scale g to be fully determined At T ^ 0, combining the OPT and RG equation fixes X{g,/T) and m{g,/T) 
for a given g/T input, which also fixes P/Po as in (7.16). But this combined solution is somewhat academic for the 
(j)^ model, specially in the symmetric phase studied here, where there is no particular physical input temperature to 
compare with. Moreover, since the scale invariance of the RGOPT pressure is still imperfect at two-loop order, the 
remnant scale dependence implies different P/Pq values for different fx/T input choices. For instance, solving similarly 
the OPT and RG equations for t = g/{2TrT) = 1/2 (respectively t = 2) gives P/Pq = 0.881 (respectively 0.921). This 
is consistent with the previously examined remnant scale dependence of the RG and OPT pressures at such relatively 
large coupling values. 

As above mentioned one may expect that the full RGOPT solution for arbitrary scale input could give a sensible 
approximation of the genuine nonperturbative scale dependence of the coupling. What Eq. (7.16) indicates is that, 
for the physically reasonable scale choice /i ~ 2ttT, the optimized coupling g and mass rh/T are both of order 1, lying 
in the rather nonperturbative range where the soft modes ~ gT become comparable with the hard modes ~ T . Had 
we rather found optimized values i, g <C 1, we would not learn much beyond standard perturbation theory. However, 
we cannot easily follow the combined solution over an arbitrarily large range of scale: the physical branch solution 
remains real for relatively large variations of t > 1, but becomes complex for t < 0.27, which corresponds to the 
already mentioned complex RG solntion for the large coupling A = —35o/(25i). Thus as already stressed we should 
not trust our results for g above 5 > 1, which corresponds to t < .5. For t > I (respectively t < 1) within a reasonable 
range, the combined RGOPT solution gives slightly smaller (respectively larger) optimized g values as compared 
with the one in (7.16), unlike the standard perturbative RG behavior of the (j)'^ model at small coupling. Indeed for 
even larger /x ^ 2ttT^ where one expects to recover the four-dimensional T = 0 model properties, the RGOPT 
real solution gives a slowly (logarithmically) decreasing coupling and mass, which appears roughly consistent naively 
with the expected triviality property[47]. Gonversely, for g <C 2'kT one expects to recover the three-dimensional 
high-T limit, such that after eventually reaching a maximum, the optimized coupling is expected to decrease again for 
g <C 2ttT. But it is difficult to follow the full RGOPT solution becoming complex for g <C 27rT. Nevertheless, even if 
they have no real intersections, the OPT and RG solutions can be reliably determined for g <C 2'kT, as discussed in 
Sec. VH.B, and for fixed m <C 2'itT, \o¥t{p/T) or X^oig/T) decrease logarithmically for g <C 2ttT. 

Note that when incorporating the S 2 0 term from Eq. (3.7), the full RGOPT solutions similar to (7.16) are shifted 
to a different scale but with the same qualitative behavior, so the net effect of S 2 0 appears essentially as a 
renormalization scheme redefinition, without drastically changing the results. Incidentally, for S 2 ^ 0 the full RGOPT 
solution is already complex for t = 1, while real solutions appear for slightly larger t > 1.5 values, with corresponding 
optimal P/Pq values very close to the one in (7.16). This simply reflects that the occurrence of complex RGOPT 
solutions in a given theory depend much on the renormalization scheme, so that real solutions could be recovered in 
principle from appropriate perturbative scheme changes [17]. But this is a much more involved task in the present 
finite temperature case. Moreover, since the non-trivial RG and OPT intersecting solution happens first at two-loop 
order, it is probably safer not to take as a very firm prediction the result in (7.16), nor to vary in a wide range around 
the preferred value g ^ 2'kT. From the example of T = 0 results in other models [17] the RGOPT variational fixed 
point solution is more likely to stabilize at the three-loop order, and probably with a more perturbative value of the 
optimized coupling. 


VIII. CONCLUSIONS AND PROSPECTS 

Let us summarize our main results. We have shown that the standard treatment of the free energy (pressure) in 
massive thermal theories, with minimally subtracted counterterms, as primarily done in resummation approaches 
like OPT/SPT,HTLpt typically, lacks RG invariance. We have then recalled a general simple recipe to restore RG 
invariance, leading unavoidably to additional finite, temperature-independent vacunm energy contributions, sys¬ 
tematically derivable in a perturbative fashion. We have next explained that the OPT/SPT,HTLpt resummation 
methods based on the modification of the perturbative expansion from the linear ^-expansion in general do not 
preserve RG invariance, which can however be restored for a different interpolating prescription, Eq. (4.1), uniquely 
dictated by universal first order RG coefficients, a = 7o/&o- Moreover, the resulting RGOPT resummation can use 


Similarly for QCD at zero temperature[16, 17] the pion decay constant has been used as a reference physical scale, and the analog of 
the combined RG and OPT solutions completely fixed F-jt/^qcd- 
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the RG equation as an alternative or additional combined constraint to determine a self-consistent nonperturbative 
mass (and coupling), in addition to the sole standard OPT prescription using only the mass optimization. We 
have then illustrated the RGOPT in details by evaluating the free energy of a thermal scalar field theory at one- 
and two-loop level. The results show a substantial improvement regarding this type of nonperturbative approach. 
Namely, we obtain exact RG/scale invariance at the first non trivial one-loop RGOPT order, which also reproduces 
all the exact large-iV results [22] of the 0{N) scalar model. At two-loop RGOPT order, the scale dependence and 
stability are drastically improved up to relatively large coupling values, with respect to most of the other available 
resummation approaches. Not surprisingly the pressure obtained from the RG mass gap equation happens to have 
better convergence and scale dependence properties for large coupling than the pressure obtained from the OPT mass 
gap. We have also illustrated the full RGOPT result obtained from combining the RG and OPT equations, therefore 
completely fixing the coupling and mass for a given input scale, with results that can be considered as a variational 
approximation to the truly scale invariant nonperturbative all order result. Beyond two-loop order, since all relevant 
perturbative results are available at three-loop order [7], the RGOPT procedure can be applied, but we leave this 
for future work as it becomes algebraically and numerically somewhat more involved. Besides, we are confident 
that it will further improve the scale dependence with respect to the two-loop results, without explicit three-loop 
calculations: as explained the RGOPT construction will guarantee the RGOPT free energy to be perturbatively RG 
invariant up to neglected four-loop terms, for arbitrary mass m. This implies, after using the (OPT or 

RG) gap equation, that perturbatively a remnant scale dependence appears in the free energy at order ~ A"*. The 
same line of reasoning also explains why the lack of RG invariance of OPT/SPT at formally one-loop order for 
arbitrary m, remains partly hidden at one- and two-loop but resurfaces maximally at perturbative three-loop order 
A^. We also see that the RGOPT anticipates the predictions by one perturbative order, the exact one-loop results 
being qualitatively similar to standard two-loop resummation results. Therefore, one may argue similarly that the 
two-loop RGOPT results should be qualitatively comparable to (standard) SPT three-loop results (with a sensibly 
better scale invariance). Indeed, incorporating the $2 subtraction term at two-loop order, which includes a (RG) 
part of the three-loop contributions, further improves the convergence and scale dependence in the nonperturbative 
coupling range. For these reasons considering the full three-loop contributions is certainly a welcome refinement but 
not necessary a crucial one in order to demonstrate the efficiency of the method which constitutes our main goal here. 

It should be also clear that the whole construction illustrated in particular in Secs. 3 to 5 is actually more general, 
and that it is applicable to QCD. We have already mentioned some properties anticipated to be similar, or sometimes 
different, in the QCD case. One could expect that the thermal QCD application may be a priori much more difficult 
than the traditionally simpler scalar model. But given that the difficult gauge dependence and related QCD issues 
have been solved by the HTL formalism [4, 8], the elaborate perturbative HTLpt calculations performed for thermal 
QCD up to three-loop order in [9] should be readily adaptable to our RGOPT method, which in a first stage relies 
entirely on perturbative calculations. In HTLpt only vacuum energy, mass, and coupling counterterms are necessary 
to renormalize the thermodynamic potential. The quark mass anomalous dimension is just the standard one, while 
the gluon mass anomalous dimension is easily extracted from corresponding known counterterms, given e.g. in [8, 9]. 
Thus our subtraction procedure to recover RG invariance will work just like in the scalar model, applying RG with 
Eq. (3.1), (3.3), and the modifed interpolation Eq. (4.1). Moreover, the HTLpt formalism is inherently a high-T 
expansion, therefore it will give OPT and RG equations in m/T and g, simpler than the somewhat involved two-loop 
exact T-dependent mass gap Eqs. (7.2), (7.3) (except that for QCD, HTLpt at two-loop order involves mlnm terms, 
due to the two-loop QCD free energy graph structure). Finally, for QCD the known first four coefficients of the 
/3-function, bo - • - b^, all have the same (negative) sign, such that no fake perturbative fixed points, such as the one 
present here in the scalar model, will obstruct the identification of correct RGOPT solutions. We are confident that 
a similarly improved scale-dependence and overall stability will be obtained also from RGOPT adaptation of HTLpt, 
which could potentially put confrontation with thermal QCD lattice results on even firmer grounds. 
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Of course with the limitations that our two-loop results, while resumming some RG dependence of higher order terms, does not reproduce 
the full three-loop contributions, in particular those given by the three-loop thermal “basketball” graph [7, 48]. When the subtraction 
term S 2 X from (3.7), (3.8) is included, it is related from RG properties to the single logarithm coefficient of this three-loop basketball 
graph. 
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Appendix A: Alternative derivation of subtraction terms from bare RG-invariance 


For completeness is this appendix we derive the subtraction terms in Eq. (3.8) alternatively by RG invariance 
considerations solely on the bare expression. The bare free energy at two-loop level consists of Eq.(2.4) plus the 
remnant divergent terms [7, 34] (after mass and coupling renormalization at this order) with Z? = 4 — 2e: 


(47r)^£lo (residual) = —m^ 


1 

4 ^ 


1 



(Al) 


As amply explained, minimally subtracting Eq.(Al) would not produce a finite RG-invariant expression. Instead, an 
explicitly RG-invariant counterterm can be written [17, 29, 30] in the general (perturbative) form in terms of the RG 
invariant bare mass and couplings as: 


m. 


An 


(A2) 


with HqY) = arbitrary series in e, most conveniently determined perturbatively order by order. Now 

demanding Eq. (A2) to cancel the remnant divergent terms in (Al), upon using in Eq. (A2) the well-known expressions 
of the mass and coupling counterterms up to two-loop order, reading in our convention: 


= 


A, 


A 



(A3) 


m 2e 


7o(7o + bo) 
8e2 


7i 

4e 


A2 + . 


(A4) 


expanded in A and e series, after some algebra it uniquely determines ho = sq and hi = si, given in Eq. (3.8). But 
it leave additional finite subtraction terms identical to Eq. (3.8). Note that it necessarily involves an rrYso/X: using 
in Eq.(A2) simply mo cannot cancel the one-loop divergence in Eq.(Al) since the latter is (!1(A°). Of course for the 
present (f)'^ model this construction is equivalent to the one performed in [34]. 
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